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SUMMARY 


A  further  examination  of  previously  reported  unsteady  aerodynamic  data 
on  a  tunnel  spanning  wing  (both  swept  and  unswept),  obtained  in  the  United 
Technologies  Research  Center  (UTRC)  Main  Wind  Tunnel,  was  performed.  In  the 
original  study  (NASA  CR-3092)  the  analysis  of  these  data  was  focused  on  the 
behavior  of  the  integrated  force  components  of  the  NACA  0012  airfoil  section 
tested  (lift,  drag,  and  pitching  moment).  The  objective  of  the  present 
investigation  was  to  study  the  unsteady  chordwise  force  response  on  the  air¬ 
foil  surface  and  to  examine  its  sensitivity  to  the  various  system  para¬ 
meters  . 

The  main  body  of  this  data  analysis  was  carried  out  by  analyzing  the 
propagation  speed  of  pressure  disturbances  along  the  chord,  and  by  studying 
the  behavior  of  the  unsteady  part  of  the  chordwise  pressure  distribution  at 
various  points  of  the  airfoil  pitching  cycle.  It  was  found  that  Mach  number 
effects  dominate  the  approach  to  stall  and  the  inception  of  stall,  both 
static  and  dynamic;  the  stall  angle  decreases  as  the  Mach  number  increases. 

It  was  also  found  that  sweep  dominates  the  load  behavior  within  the  stall 
regime;  large  phase  differences  between  unswept  and  swept  responses,  that  do 
not  exist  at  low  lift  coefficient,  appear  once  the  stall  boundary  has  been 
penetrated. 

Conventional  pressure  normalization  procedures  were  reviewed  and  it  was 
found  that  the  cosine  law  for  sweep  angle  is  not  applicable  at  angles  of 
attack  above  the  stall  angle,  both  for  steady  and  unsteady  conditions.  It 
was  also  found  that  the  reduced  frequency  parameter  traditionally  used  in 
unsteady  aerodynamic  analysis  is  not  a  reliable  indicator  of  the  unsteady 
aerodynamic  response  in  the  high  angle  of  attack  regime. 

An  additional  analysis  was  performed  to  resolve  a  troublesome  aspect  of 
the  previously  reported  results.  It  was  found  that  the  relative  displacement 
of  the  lift  and  moment  hysteresis  loop  response  (excluding  nonoscillatory 
terms)  during  upstroke  between  unswept  and  swept  wings  is  a  real  effect,  and 
is  almost  entirely  described  by  the  phase  and  amplitude  differences  in  the 
first  harmonic  components. 


INTRODUCTION 


The  need  for  swept  wing  unsteady  aerodynamic  data  into  the  dynamic  stall 
regime  was  addressed  in  Ref.  1  wherein  test  results  were  reported  from  exper¬ 
iments  performed  using  a  NACA  0012  airfoil  section  tunnel-spanning  wing  (TSW) 
oscillating  in  pitch  in  both  the  unswept  and  30  deg  swept  back  configura¬ 
tions.  The  main  motivation  for  obtaining  such  data  is  for  use  by  the 
helicopter  designer  who  must  account  for  a  widely  varying  aerodynamic 
environment  on  the  rotor  blade  as  it  traverses  the  azimuth  plane.  In  addi¬ 
tion  to  this  immediate  engineering  application,  these  results  provide  the 
opportunity  to  examine  the  effects  of  sweep  at  moderate  pitching  amplitudes, 
and  thus  permit  the  evaluation  of  techniques  currently  in  use  for  estimating 
the  unsteady  swept  wing  aerodynamic  forces  from  the  existing  data  banks  for 
zero  sweep  angle. 

The  results  that  were  reported  in  Ref.  1  focused  on  (1)  the  unsteady 
integrated  load  components  (lift,  drag,  and  moment)  computed  from  unsteady 
chordwise  pressure  distributions  for  both  stalled  and  unstalled  conditions, 
and  (2)  the  surface  pressure  fluctuations  obtained  in  the  stall  regime  from 
companion  steady  airfoil  experiments.  Several  findings  were  reported,  of 
which  the  most  notable  is  that  sweep  tends  to  delay  the  onset  of  dynamic 
stall.  A  summary  of  the  main  results  can  be  found  on  page  40  therein.  A 
companion  data  report  to  Ref.  1  was  also  published  (Ref.  2)  containing  a 
complete  set  of  unsteady  normal  force,  chord  force,  lift  force,  drag  force, 
and  moment  responses  in  the  form  of  hysteresis  loops.  The  data  report  was 
arranged  to  emphasize  the  effect  of  sweep  on  the  results.  A  review  of 
related  work  in  this  field  can  be  found  in  the  introduction  to  Ref.  1  and  in 
the  associated  bibliography. 

The  present  report  utilizes  the  same  data  base  that  was  used  to  produce 
the  results  in  Refs.  1  and  2  and  is  essentially  a  continuation  of  that  data 
analysis.  Here,  the  emphasis  is  on  the  unsteady  chordwise  pressure  data  and 
their  sensitivity  to  the  various  test  parameters:  wing  sweep  angle,  A; 
pitching  amplitude,  a;  mean  angle  of  attack,  aM;  pitching  frequency,  f;  and 
Mach  number,  Mc  *  cos  A.  Two  parametric  studies  are  presented  herein: 

(1)  the  behavior  of  the  chordwise  propagation  of  pressure  waves  associated 
with  the  dynamic  stall  process,  and  (2)  the  behavior  of  the  unsteady  surface 
pressures,  both  individually  as  time  histories  over  the  whole  motion  cycle, 
and  collectively  as  chordwise  distributions  at  specific  points  of  the  motion 
cycle. 

An  accompanying  data  report  (Ref.  3)  contains  listings  of  the  first  ten 
harmonic  components  of  the  individual  chordwise  pressure  responses,  and  of 
the  integrated  load  results.  The  data  tables  in  Ref.  3  can  be  used  to  repro¬ 
duce  the  oscillatory  components  of  the  individual  pressure  time  histories, 
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and  the  complete  normal  force,  chord  force,  lift  force,  drag  force,  and 
moment  coefficient  time  histories.  (Here,  the  term  "complete”  means 
nonoscillatory  plus  oscillatory  parts,  and  it  is  understood  that  the  chord 
and  drag  forces  were  computed  by  integration  of  pressure  distributions  and 
therefore  do  not  include  the  effects  of  viscous  drag.)  Hysteresis  plots  of 
the  ingegrated  force  coefficient  results  (versus  angle  of  attack)  were 
reported  in  Ref.  2. 

The  present  document  begins  with  a  review  of  the  data  base  and  the 
approach  to  the  normalization  of  the  data.  This  is  followed  by  a  description 
of  the  data  retrieval  process  that  was  used  to  expedite  the  performance  of 
the  parametric  analysis.  A  discussion  of  the  chordwise  wave  propagation 
phenomenon  and  its  sensitivity  to  the  test  parameters  is  then  presented.  A 
review  of  the  effectiveness  of  the  data  normalization  procedure  is  presented 
next.  The  report  then  proceeds  with  the  analysis  of  the  chordwise  response 
behavior  to  the  various  test  parameters,  including  a  brief  examination  of  the 
individual  pressure  time  histories  very  near  the  leading  edge.  Finally,  a 
troublesome  aspect  of  the  hysteresis  loop  comparisons  of  Ref.  1  is  revisited 
and  resolved. 
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LIST  OF  SYMBOLS 


cp(x,t) 


p(x  ,t) 


dimensionless  pitch  rate,  Eq.  (14) 

dimensionless  pitch  rate  at  vortex  inception  angle 

sound  speed 

airfoil  chord 

chord  force  coefficient 

drag  force  coefficient,  Eq.  (5) 

lift  force  coefficient,  Eq.  (4) 

moment  coefficient 

normal  force  coefficient 

pressure  coefficient  distribution,  Eq.  (3) 
airfoil  pitching  frequency,  Hz 
chordwise  reduced  frequency,  Eq.  (2) 
chordwise  Mach  number,  Eq.  (1) 
free  stream  Mach  number 
slope,  Eq.  (11) 

chordwise  pressure  distribution,  newton/m 
dynamic  pressure  based  on  chordwise  velocity,  newton/m 
period,  sec 
time,  sec 

vortex  inception  time,  sec 
chordwise  velocity,  ra/sec 
free  stream  velocity,  m/sec 


4 


LIST  OF  SYMBOLS  (Cont'd) 


V 


W 


a 

a 

a 

a 

a 

a 

e 


M 

s 

ss 

V 

vo 


A 


T 


T 


O 


X 

(0 


average  dimensionless  wave  speed,  Eq.  (8) 

average  dimensionless  wave  speed  between  stations  i  and  j 
airfoil  angle  of  attack,  deg,  Eq.  (6) 
mean  angle  of  attack,  deg 

dynamic  stall  inception  angle,  deg,  Eq.  (19) 

steady  state  stall  angle,  deg,  Eq.  (19) 

vortex  inception  angle,  Eq.  (10) 

constant  used  in  Eq.  (11) 

data  scatter  parameter,  Eq.  (15) 

sweep  angle,  deg 

dimensionless  time  based  on  oscillation  period,  also 
fraction  of  period 

dimensionless  vortex  inception  time  based  on  oscillation  period 

phase  lead,  deg 

distance  from  airfoil  nose 

circular  frequency,  rad/sec 
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DESCRIPTION  OF  THE  DATA  BASE 


The  data  base  for  the  present  study  was  obtained  from  aerodynamic  exper¬ 
iments  performed  on  an  oscillating  tunnel-spanning  wing  (TSW)  in  both  unswept 
and  30  deg  swept  configurations  in  the  United  Technologies  Research  Center 
(UTRC)  Main  Wind  Tunnel.  The  chordwise  profile  taken  normal  to  the  wing 
leading  edge  is  a  NACA  0012  airfoil.  A  detailed  description  of  the  test 
program  can  be  found  in  Ref.  1.  Although  the  data  base  and  normalization 
procedure  are  also  fully  described  therein,  the  portions  that  are  pertinent 
to  the  present  study  are  repeated  here  for  completeness.  Of  primary  interest 
in  this  investigation  is  the  unsteady  chordwise  pressure  distribution.  The 
unsteady  data  acquisition  points  that  were  used  along  the  chord  are  schem¬ 
atically  depicted  in  Fig.  1. 

The  experiments  of  Ref.  1  were  performed  for  the  matrix  of  test  para¬ 
meter  settings  shown  in  Table  1.  These  parameters  were  the  wing  sweep  angle, 
A,  the  amplitude  of  pitching  motion,  a,  the  mean  angle  of  attack,  a^,  the 
pitching  frequency,  f,  and  the  approach  Mach  number,  ^ .  The  effect  of 
sweep  was  accounted  for  in  the  normalization  process  by  referring  all  of  the 
data  to  the  component  of  Mach  number,  Mc,  normal  to  the  wing  span.  This 
quantity  is  related  to  the  freestream  value,  ^  ,  by  the  formula 

Mc  =  cos  A  ( l ) 

where  subscript  c  denotes  the  chordwise  direction.  Data  were  obtained  for 
two  values  of  Mc  (=  0.30  and  0.40)  for  each  value  of  A  (=0  deg  and  30 
deg).  This  approach  required  testing  the  30  deg  swept  configuration  at  a 
higher  freestream  velocity  as  shown  in  Table  1.  Unsteady  pressure  data  were 
obtained  for  two  values  of  a  (8  deg  and  10  deg)  at  four  mean  angle  of  attack 
settings  (0,  9,  12  and  15  deg)  about  the  quarterchord  axis.  The  angles  were 
measured  in  a  plane  perpendicular  to  the  span  axis.  The  choice  of  pitching 
frequencies  used  during  the  test  program  was  predetermined  according  to  Table 
2  which  lists  the  nominal  values  of  reduced  frequency  defined  by 


k 


c 


TT  cf 


(2) 


where  is  the  speed  of  sound. 
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The  effects  of  sweep  on  the  pressure  data  were  accounted  for  in  Ref.  1 
by  normalizing  the  data  with  respect  to  the  velocity  component  normal  to  the 
span  to  obtain  the  following  expression  for  dimensionless  pressure 

Cp  (x,t)  =  -  p(x ,t)/qcos2A  (3) 

where  x  is  the  dimensionless  position  along  the  chord  and  t  is  time.  This 
nondimensionalization  scheme  is  based  on  a  steady-state  sweep  analysis  which 
assumes  the  flow  to  be  potential  and  the  pressure  distribution  on  the  wing  to 
be  determined  entirely  by  the  magnitude  of  the  normal  velocity  component, 

V^,  cos  A.  Note  also  that  in  this  report,  the  pressure  coefficient  is 
defined  as  the  negative  of  the  normalized  pressure.  This  is  convenient  for 
relating  positive  integrated  lift  to  a  positive  value  of  Cp. 
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DATA  PROCESSING  AND  ANALYSIS 


Introductory  Comments 

The  analysis  of  unsteady  chordwise  pressures  and  their  sensitivity  to 
the  various  system  parameters  by  conventional  methods  is  a  formidable  under¬ 
taking  for  a  data  base  the  size  of  the  one  shown  in  Table  1.  This  is  true 
mainly  because  an  important  aspect  in  the  study  of  chordwise  pressure 
behavior  involves  a  comparative  study  of  the  pressure  distribution  at  various 
points  of  the  airfoil  pitching  motion  cycle.  This  approach  requires  substan¬ 
tial  data  manipulation  if  a  reasonable  number  of  representative  cases  and 
comparisons  are  to  be  analyzed  in  any  significant  detail.  (It  is  noted  that 
the  chordwise  pressure  in  the  present  study  consists  of  21  channels  of  data, 
including  13  upper  surface  and  8  lower  surface  locations.) 

The  approach  that  was  used  to  expedite  the  present  analysis  consisted  of 
(1)  compacting  the  processed  data  into  a  data  file  where  each  record  (test 
case)  is  addressable  either  by  record  number  or  by  the  corresponding  test 
parameters  (A,  a,  a^,  f,  and  Mc)  that  are  unique  to  the  record,  and  (2) 
tailoring  a  data  management  program  where  specific  channels  from  any  combina¬ 
tion  of  records  can  be  retrieved  for  numerical  and/or  graphical  analysis  by 
means  of  a  series  of  built-in  computer  queries. 


Data  Compaction  and  Storage 

The  direct  access  data  file  that  was  used  in  the  present  analysis  was 
generated  in  two  steps:  (1)  compacting  the  unswept  and  swept  data  bases  into 
a  common  data  file,  and  (2)  creating  a  file  directory  for  interfacing  the 
data  file  with  the  data  retrieval  software  (described  later  in  this 
sect  ion) . 

Data  compaction  was  achieved  for  the  unswept  and  swept  data  bases  by 
cycle-averaging  the  individual  pressure  responses  of  each  case  over  5  cycles 
of  the  airfoil  pitching  motion  and  then  reducing  the  resulting  time  histories 
to  24  harmonic  components  (including  the  zeroth  harmonic  term).  The 
fundamental  frequency  used  in  the  Fourier  analysis  process  was  the  pitching 
frequency  of  the  airfoil.  The  unswopt  and  swept  wing  data  bases  were  then 
combined  into  a  common  compacted  data  file. 

The  normal  force,  chord  force,  and  moment  responses  were  then  directly 
computed  from  separate  chordwise  integrations  (by  trapezoidal  rule)  of  the 
in-phase  and  out-of-phase  portions  (relative  to  the  airfoil  motion)  of  each 
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harmonic  component.  The  lift  and  drag  forces  were  then  computed  from  the 
normal  force  and  chord  force  using  the  following  formulas: 


CL  =  CN  cos  0  +  Cc  s*n  a  (A) 


CD  =  CN  s*n  a  “  Cc  cos  a  (5) 

where  a  is  given  by 

a  =  +  a  sin  «t.  (6) 

A  schematic  depicting  the  force  vectors  is  shown  in  Fig.  2.  Bv  inspection  of 
these  equations,  it  is  seen  that  the  harmonic  coefficients  of  and  CD 

cannot  be  computed  directly  from  the  harmonic  coefficients  of  and  Cc 

because  of  the  nonlinear  appearance  of  a  within  the  trigonometric  terras. 
Therefore,  an  intermediate  computatin  of  the  and  responses  was  made 
by  substituting  time  history  reproductions  of  C^,  Cc,  and  a  into  Eqs.  (A) 
through  (6).  The  resulting  and  Cq  time  histories  were  then  Fourier 
analyzed,  and  together  with  C^,  Cc,  and  C^,  were  added  to  the  compacted 
data  file  by  extending  the  number  of  channels  of  each  test  case  from  22  (a 
and  21  chordwise  pressures)  to  27. 

The  final  step  in  completing  the  contents  of  the  compacted  file  was 
adding  to  each  case  the  steady  state  force  and  moment  coefficients 
corresponding  to  CN,  Cc,  CM,  C^,  and  CD  as  well  as  the  actual  test  parameter 
values  (A,  a,  aM,  f,  and  Mc)  measured  during  the  test. 

The  interface  of  the  data  file  with  data  retrieval  software  was  achieved 
by  installing  a  file  directory  header  into  the  data  file.  The  file  directory 
was  indexed  to  match  the  actual  record  sequence  of  the  file  and  was  set  up 
with  nominal  values  of  A>  a,  aM,  f,  and  Mc  corresponding  to  the  actual  test 
parameter  values  of  each  record.  The  objective  here  is  that  all  requests  to 
the  data  file  are  processed  by  means  of  the  file  directory.  Data  retrieval 
is  thus  achieved  without  having  to  directly  search  through  the  individual 
records  of  the  file.  This  is  an  especially  desirable  approach  when  records 
are  being  retrieved  out  of  sequence  from  the  file. 


Data  Analysis  Techniques 

The  objective  of  the  data  analysis  in  this  report  is  to  examine  and 
evaluate  the  effect  of  the  test  parameters  (A,  a,  a^,  f,  and  Mc)  on  the 
unsteady  pressure  response  over  the  surface  of  the  blade.  To  accomplish 
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this,  the  authors  found  it  necessary  to  use  several  graphical  representations 
while  concentrating  on  the  various  characteristic  features  of  the  data.  It 
is  appropriate,  therefore,  that  the  descriptions  of  each  figure  type  be 
consolidated  in  the  present  section  so  that  subsequent  sections  can  be 
devoted  primarily  to  the  discussion  of  the  physical  phenomena  being 
depicted . 

Much  of  this  report  deals  with  changes  in  unsteady  blade  pressure 
associated  with  changes  in  certain  of  the  test  parameters.  Recall  that  in 
this  report,  pressure  coefficient  is  defined  as  the  negative  of  the 
normalized  pressure  at  any  station, 

Cp  (x.f)  =  -p(x»t  )/qcos2A  .  (7) 


In  some  figures  the  pressure  coefficient  is  plotted  versus  dimensionless 
chord  position,  Xi  while  in  others  the  independent  variable  is  dimensionless 
time,  t.  For  example,  Fig.  3  represents  the  pressure  coefficient  plotted 
versus  x  at  the  instant  that  0=0^  (upper  panel),  and  at  the  peak 
incidence  angle,  a  =  +  a  (lower  panel).  The  direct  comparison  to  be  made 

here  between  the  solid  and  dashed  lines  is  between  two  chordwise  Mach 
numbers.  (Other  plots  will  show  the  effect  of  other  parameters  in  the  same 
manner.)  A  similar  chordwise  pressure  plot,  shown  in  Fig.  4,  depicts  the 
first  harmonic  pressure  amplitude  (upper  panel)  and  phase  angle  (lower  panel) 
taken  from  the  Fourier  representation  of  the  pressure  time  history.  Here  the 
direct  comparison  is  between  the  two  sets  of  symbols  representing  the  two 
sweep  angles  tested. 

A  third  type  of  chordwise  pressure  plot  is  shown  in  Fig.  5.  Here,  the 
format  of  the  individual  panels  is  similar  to  that  of  Fig.  3,  but  the  intent 
of  the  figure  as  a  whole  is  significantly  different.  Now  the  solid  and 
dashed  lines  denote  upper  and  lower  surface  pressures,  and  the  major  compar¬ 
ison  is  made  between  the  two  columns.  (In  this  case,  the  effect  of  sweep  is 
depicted.)  A  second  parametric  comparison  can  be  made  in  Fig.  5  by  viewing 
successive  pairs  of  panels  which,  in  this  case,  are  for  three  instantaneous 
angles  of  attack,  a.  Both  the  numerical  value  of  a  and  the  relative  position 
of  the  blade  during  the  motion  cycle  are  shown  between  each  pair. 

The  last  type  of  chordwise  pressure  plot  shown  in  Fig.  6  differs  from 
Fig.  5  only  in  the  orientation  of  the  individual  plots.  Here,  for  example, 
the  upper  row  curves  were  obtained  at  a  mean  angle  of  attack  of  =  9  deg, 
and  the  lower  row  at  aM  **  12  deg.  Each  column  then  represents  a  specific 
angular  displacement  from  the  mean  (i.e.,  7,  7.3,  and  8  deg)  as  the  peak 
angle  of  attack  is  approached  and  attained  during  the  upstroke. 

Individual  pressure  time  histories  are  also  used  in  this  report.  A 
typical  example  is  shown  in  Fig.  7.  In  every  case  the  upper  row  of  panels 


10 


will  represent  the  motion  of  the  blade  about  the  mean  angle  of  attack,  and 
each  successive  row  will  show  the  pressure  time  history  for  the  specific 
chordwise  measurement  stations  denoted  by  x  •  As  in  Fig.  5,  the  comparisons 
in  Fig.  7  are  made  between  solid  and  dashed  lines  within  each  panel,  and 
between  adjacent  columns. 


In  Ref.  4  it  was  shown  that  plots  of  unsteady  pressure  versus  both  time 
and  chordwise  position  could  be  analyzed  to  yield  quantitative  information  on 
the  chordwise  propagation  of  pressure  waves.  That  technique  is  applied  here¬ 
in.  Plots  of  pressure  time  histories  were  produced  at  several  chordwise 
positions,  and  in  the  example  of  Fig.  8,  the  data  are  restricted  to  the  first 
15  percent  of  the  chord  for  clarity.  In  this  figure  the  vertical  axis 
represents  both  the  unsteady  pressure  coefficient  level  relative  to  steady 
state  (left  scale)  and  the  chordwise  location  of  the  pressure  measurement 
(right  scale).  The  result  is  a  pseudo  three-dimensional  plot  of  the 
Cp(x»T)  surface,  referred  to  as  a  pressure  carpet  plot.  To  interpret  these 
figures,  refer  each  curve  to  its  zero  level  at  the  appropriate  x-tic  mark  on 
the  right  scale,  and  use  the  left  scale  as  a  gage  to  measure  the  pressure 
relative  to  its  zero  “vel.  The  precipitous  drop  in  C in  the  vicinity  of 
t*=0.17  to  0.2  is  customarily  associated  with  dynamic  stall  and  is  usually 
ascribed  to  the  formation  and  chordwise  passage  of  a  stall  cell  or  vortex. 


For  this  carpet  plot  a  set  of  constant  pressure  contours  (isobars)  was 
constructed  in  the  x .T  plane  as  shown  in  Fig.  9.  The  two  columns  of  numbers 
on  the  right  represent  the  percent  of  the  pressure  range  from  minimum  to 
maximum  and  the  associated  value  of  unsteady  pressure  coefficient  relative  to 
steady  state  (for  example,  a  value  of  Cp  =  -3.0  is  4  percent  of  the  range 
above  minimum  Cp,  and  a  value  of  Cp  =  6.0  is  74  percent  of  the  range 
above  minimum).  The  numbers  that  are  appended  to  each  contour  are  the 
corresponding  values  of  Cp.  As  shown  in  Ref.  4,  the  direction  of  the 
family  of  contour  extremes  identifies  the  direction  of  wave  propagation,  and 
a  simple  estimate  of  the  slope  yields  the  wave  speed. 

The  sensitivity  of  the  contours  to  relatively  small  timewise  discrep¬ 
ancies  in  data  from  adjacent  measuring  stations  revealed  a  minor  error  in  the 
data  set,  which  is  described  in  the  Appendix.  The  error  manifested  itself  in 
small  time  or  phase  displacements  of  adjacent  signals,  and  is  occasionally 
noticeable  as  a  zigzagging  ripple  in  the  chordwise  pressure  distribution,  as, 
for  example,  in  the  lower  right  panel  of  Fig.  5,  or  in  several  panels  of 
Fig.  6. 

As  noted  in  the  Appendix  the  entire  data  set  was  uniformly  corrected  for 
the  computations  made  in  the  wave  speed  analysis,  but  the  corrections  were 
confined  to  that  section  only.  The  remainder  of  the  analysis,  including  all 
of  the  tabulated  data  in  Ref.  3,  remained  uncorrected,  to  preserve  its 
commonality  with  previous  results  from  Refs.  1  and  2. 
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Normalization  Procedure 


In  the  present  analysis  normalization  was  achieved  by  dividing  the 
negative  of  the  measured  unsteady  pressures  by  the  dynamic  pressure, 

0.5  pV^,  and  by  multiplying  the  circular  frequency  of  airfoil  pitching 
oscillation  by  c/2Vc.  Here,  Vc  is  the  chordwise  component  of  the 
freestream  velocity  (V^,  cos  A)  and  c  is  the  airfoil  chord. 
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PRESSURE  WAVE  STUDY 


Wave  Speed  Analysis 


The  use  of  pressure  coefficient  carpet  and  contour  plots  for  wave  speed 
analysis  was  introduced  in  the  previous  section,  using  Figs.  8  and  9  for  the 
forward  15  percent  portion  of  the  airfoil  as  examples.  The  correspond i ng 
carpet  and  contour  plots  for  the  entire  chord  are  presented  in  Figs.  10  and 
11.  (In  this  particular  case,  A  =  0  deg,  a  =  8  deg,  aM  =  9  deg,  f  =  10  Hz, 

Mc  =  0.3  and  kc  =  0.125.)  In  both  contour  plots  only  positive  values  of 
0^  (cf.  Eq.  7)  have  been  included,  to  focus  on  the  region  of  maximum  normal 
force  and  hence  on  the  stall  wave  propagation. 


Figure  12  is  a  contour  map  of  the  pressure  surface  in  Fig.  8,  focused  on 
the  first  four  tenths  of  the  period,  0  <  x  <  0.4.  In  this  figure,  the 
connection  between  the  contour  extremum  at  the  leading  edge  (x0  *  0.004)  with 
the  extremum  at  the  15  percent  chord  location  ( X2  =  0.149)  by  the  dashed  line 
labeled  "overall  15 Z"  is  used  for  estimating  the  wave  speed  over  the  forward 
15  percent  of  the  blade  chord.  It  is  seen  that  this  set  of  contour  extremes 
moves  toward  increasing  time  as  the  chord  is  traversed  from  front  to  rear. 

The  dimensionless  wave  speed  was  shown  in  Ref.  4  to  be  related  to  the  ratio 
of  chord  extent  to  elapsed  time  by  the  formula 


(8) 


and  with  t  =  T  T  *  2ht/oj,  and  kc  =  cu>/2Vc,  this  becomes 


V 


W 


(9) 


A  careful  study  of  the  contours  in  Fig.  12  indicates  that  the  locus  of 
contour  extremes  is  not  linear  but  can  be  represented  by  straight  line 
segments  of  different  slopes.  A  two-segment  estimate  for  this  case  is 
illustrated  in  Fig.  13.  The  end  points  of  the  first  segment  are  denoted  by 
(X0  1  T0)  and  (xi  »  t*)  and  those  of  the  second  segment  by  (xj  ,  Tj )  and  (x2  , 
t2).  Hence,  there  are  actually  two  distinct  regions  within  the  first  15 
percent  of  the  chord  that  have  different  wave  speeds.  For  this  particular 
example,  these  values  are 

(X  ,  X  ,  X  )  -  (0.004,  0.028,  0.149) 

0  1  2 

(t  ,  x  ,  t  )  —  (0.123,  0.155,  0.187) 

0  l  2 
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and  with  k  =  0.125,  application  of  the  wave  speed  Eq.  (9)  yields 


Vw  =  —•  ~  -  0.035  (first  segment) 

0  1  10 

X  -  X 

V^w  =  — =-^—  =  0.146  (second  segment) 

12  2  1 

X  -  x 

VW  mk  =  0.090  (overall,  15%) 

02  2  0 


where  the  wave  speed  subscripts  match  the  subscripts  of  the  initial  and  end¬ 
point  coordinates  of  the  segment  in  question.  Furthermore,  Fig.  14  for  the 
entire  chord  shows  that  another  discontinuity  can  be  defined  at  the  15 
percent  chord  station,  separating  the  forward  region  from  the  aft  region  of 
the  chord,  implying  additional  speed  changes  along  the  chord.  This  last 
coordinate  point  is  given  by 

(X  ,  T  )  =  (0.971,  0.284) 

3  3 

and  two  additional  wave  speeds  can  be  defined. 


0.326  (third  segment) 

0.239  (overall,  100%) 


For  convenience,  the  individual  segment  wave  speeds  for  this  example  are 
collected  below. 


Vw  =  0.035 
01 


Vw  >  0.090 
02 


V„  *  0.146  V„  *=  0.239 

12  03 

Vw  =  0.326 
23 


It  is  seen  that  the  wave  is  slowest  in  the  region  of  high  pressure  gradient 
near  the  leading  edge  (V„  ).  This  is  the  region  characterized  by  the  dense 
array  of  constant  level  contour  lines  in  Fig,  13.  Conversely,  the  wave 
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travels  more  rapidly  over  the  second  segment  (Vw  )  where  the  chordwise 
pressure  gradient  is  lower,  and  attains  its  highlit  speed  over  the  last 
segment  )  where  pressure  gradient  is  minimal.  However,  without 

additional  ^idence,  it  is  impossible  to  establish  a  cause-and-ef feet 
relationship  between  the  pressure  gradient  and  the  speed  variation. 
Alternatively,  starting  inertia  in  the  vortex  rollup  associated  with  dynamic 
stall  could  be  the  reason  for  wave  acceleration.  Whatever  the  cause,  a 
consistent  observation  has  been  made  (as  will  be  shown  presently),  that  the 
wave  speed  seems  to  be  inversely  related  to  the  unstalled  chordwise  pressure 
gradient . 

Before  proceeding  with  a  detailed  examination  of  parametric  trends, 
several  values  of  wave  speed  for  the  unswept  wing  tests  are  compared  in  Fig. 
15  with  those  obtained  in  Ref.  4.  Overall  averaged  wave  speed  for  the  entire 
chord  (l?w  ),  and  local  averaged  wave  speed  for  15  percent  of  the  chord 

(Vu  ),  a^e  plotted  versus  reduced  frequency  in  the  upper  and  lower  panels  of 

”0  2 

Fig.  15,  respectively.  The  original  data  are  denoted  by  the  open  symbols, 
and  the  current  data  for  a  =  8  deg  and  Mc  =  0.3  by  the  solid  symbols. 

It  is  seen  that  the  current  results  are  consistent  with  previous  calcula- 
t ions . 

A  study  was  performed  to  examine  the  effect  of  each  of  the  test  para¬ 
meters  on  wave  speed  and  the  results  are  shown  in  Figs.  16  through  18.  In 
each  figure  the  overall  wave  speed,  ,  is  plotted  versus  reduced  fre- 

quency,  kc .  Mean  angle  of  attack  variation  is  always  depicted  by  the  three 
symbols  used:  circle,  triangle,  and  square.  In  each  figure  there  are  four 
panels,  with  one  parameter  held  constant  from  left  panel  to  right  panel,  a 
second  parameter  held  constant  from  top  panel  to  bottom  panel,  and  the  varia¬ 
tion  of  the  parameter  of  primary  interest  for  each  figure  highlighted  by  the 
contrast  between  open  and  solid  symbols.  In  Fig.  16  the  effect  of  amplitude 
variation  is  examined  by  noting  that  the  open  symbols  represent  a  =  8  deg 
and  the  solid  symbols  represent  a  =  10  deg.  This  figure  shows  little  or 
no  effect  of  amplitude  on  overall  wave  speed.  Although  it  was  originally 
suggested  in  Fig.  20  of  Ref.  4  that  increasing  amplitude  appeared  to  increase 

the  value  of  Vu  at  any  k  ,  the  present  results  do  not  corroborate  this 
"0  3  .  c 

earlier  indication.  It  is  noteworthy  that  the  data  in  Ref.  4  were  obtained 
for  a  variety  of  Mach  number  values  and  of  profile  shapes,  whereas  the  pre¬ 
sent  results  are  taken  from  a  wholly  self-consistent  data  base.  Also,  this 
figure  and  subsequent  figures  show  that  there  is  little  or  no  effect  of  mean 
angle  of  attack  on  wave  speed,  which  agrees  with  the  observations  made  in 
Ref.  4. 


Mach  number  is  the  parameter  of  primary  interest  in  Fig.  17,  and  its 
effect,  although  small,  seems  to  be  consistent.  For  A  ■  0  deg  (top  panels) 
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the  higher  Mach  number  (solid  symbols)  yields  a  slightly  smaller  value  of 
Vy  than  Mc  =  0.3,  although  the  two  sets  of  data  appear  to  converge  as  k 
increases.  For  A  =  30  deg  (bottom  panels)  there  is  less  distinction  between 
the  two  sets  of  Mach  number  data,  and  it  could  be  argued  that  the  effect  of 
Mach  number  is  minimal. 


By  far,  the  parameter  that  has  the  greatest  effect  on  Vy  is  the  sweep 
angle,  as  shown  in  Fig.  18.  In  all  four  panels  the  overall  wave  speed  for 
the  unswept  wing  (open  symbols)  is  consistently  and  significantly  greater 
than  that  for  the  swept  wing  (solid  symbols)  by  a  factor  of  1.5  to  2.0.  This 
phenomenon  of  lower  chordwise  wave  speed  in  swept  flow  has  not  been  predicted 
by  any  previous  analysis  nor  observed  in  any  previous  experiment. 
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In  an  effort  to  examine  these  trends  more  closely,  wave  speeds  for  the 
several  segments  (defined  earlier)  were  computed  and  are  displayed  in  Figs. 
19  and  20  for  Mc  =  0.3  and  0.4,  respectively  (at  the  same  value  a  =  8 
deg).  The  three  left  hand  panels  in  each  figure  are  the  three  local  wave 
speeds,  V. 


01 


*W, 


and  V, 


,  proceeding  aft  along  the  chord.  The  two  right 
'  w23  _ 

hand  panels^are  the  two  average  wave  speeds,  Vu  over  the  first  15  percent 

of  the  chord,  and  Vy  over  the  entire  chord.  (The  latter  are  repeated  in 

these  two  figures  from  the  top  panels  of  Fig.  18.)  It  is  clearly  shown  here 

that  there  is  substantially  no  effect  of  sweep  angle  on  the  local  wave  speed 

over  the  forward  15  percent  of  the  chord,  either  by  examination  of  the  panels 


for 


and  V, 


Wq i  W 1 2  '  .  ~  WQ2  " 

sweep  effect  is  perceived  to  occur  in  the  aft  portion  of  the  chord  (lower 

left  panel  of  both  Figs.  19  and  20),  with  factors  of  as  much  as  2  observed  in 

swept  wing  wave  speed  versus  that  for  the  unswept  wing.  This  result  is 

representative  of  all  cases  examined.  If  the  pressure  wave  is  associated 

with  the  formation  and  chordwise  transport  of  a  vortex,  these  results  would 

imply  that  the  vortex  motion  along  the  chord  is  independent  of  sweep  over  the 

first  15  percent  region  of  the  chord,  and  is  strongly  dependent  on  sweep 

thereafter . 


or  from  the  results  for  V, 


Wn 


Virtually  all  of  the 


Finally,  an  analysis  of  the  chordwise  variation  of  local  wave  speed  is 
provided  by  Fig.  21  for  the  case  with  Mc  =  0.4,  a  =  8  deg  and  =  15  deg. 

In  this  figure  the  local  wave  speed  is  plotted  versus  dimensionless  chord 
position,  Xi  for  several  values  of  k.  The  unswept  results  are  plotted  in  the 
left  panels  and  the  swept  results  in  the  right.  Linear  scales  are  used  in 
the  upper  panels,  and  two  effects  clearly  emerge.  The  first  is  the  general 
increase  of  Vy  with  k  at  any  given  value  of  x  (although  there  is  some 
scatter)  for  both  values  of  A.  The  second  is  the  continual  increase  with  x 
of  Vy  for  A  «  0  deg,  and  the  much  smaller  increase  with  x  of  Vy  aft  of  x  = 

0.1  for  A  ■  30  deg.  In  an  effort  to  spread  out  the  many  points  plotted  at 
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small  values  of  x>  the  same  data  are  repeated  in  the  semi-log  plots  in  the 
lower  two  panels.  While  there  is  no  unusual  revelation  afforded  by  these 
semi-log  plots,  they  do  appear  to  resolve  the  cluttered  points  into  an 
orderly  chordwise  progression. 


Vortex  Inception  Angle  Variation 

In  pursuing  further  the  interpretation  of  the  phenomenon  just  described 
as  the  formation  of  a  vortex  and  its  subsequent  chordwise  transport,  the 
inception  time  can  be  studied  and  related  to  the  vortex  inception  angle, 
ay.  A  direct  computation  of  ay,  can  be  made  by  substituting  the  dimen¬ 
sionless  inception  time,  xQ  =  to t Q / 2n  (cf.  Fig.  13  or  14),  into  Eq.  (6) 
for  instantaneous  angular  position, 

av  =  ajj  +  a  sin  2n0  (10) 

The  results  are  shown  in  Figs.  22  and  23  as  plots  of  vortex  inception  angle 

versus  reduced  frequency.  As  before,  the  use  of  circles,  triangles,  and 
squares  denote  aM  =  9,  12,  and  15  deg,  and  open  and  solid  symbols  highlight 
the  comparison  of  sweep  angle  effects  in  Fig.  22  and  of  Mach  number  effects 

in  Fig.  23.  These  plots  also  contain  information  on  the  steady  state  stall 

angle  (triangular  symbols  pointing  to  the  appropriate  values  on  the  ordinate) 
and  on  the  best  linear  fit  to  the  data.  The  latter  were  obtained  by  a  least 
squares  procedure  for  each  data  group,  and  the  former  by  estimating  the 
steady  state  angle  of  attack  at  which  dCM/da  *  0  in  Figs.  14  through  17  in 
Ref.  1. 


In  both  figures  the  vortex  inception  angle,  av,  varies  linearly  with 
kc,  as  indicated  by  the  good  correlation  of  the  linear  least  squares  line 


a 


v 


+  mkc 


(11) 


with  the  measured  data,  and  the  narrow  scatter  range  for  each  data  set.  (In 
this  equation  aVo  is  the  zero  frequency  intercept  and  m  is  the  slope.) 

Figure  22  shows  that  sweep  has  a  tendency  to  delay  vortex  inception,  but  this 
delay  is  less  than  two  degrees  overall,  and  less  than  one  degree  in  the 
majority  of  cases  examined.  In  contrast  to  this  result,  Fig.  23  shows  a 
substantial  Mach  number  effect,  with  incipient  vortex  angle  differences  of  as 
much  as  three  or  four  degrees  as  Mc  varies  from  0.3  to  0.4.  Both  figures 
indicate  that  the  steady  state  stall  angle  based  on  zero  moment  slope, 
ag8,  is  unaffected  by  sweep  angle  changes.  (However,  Figs.  14  through  17 
of  Ref.  1  show  that  the  shape  of  the  lift  curve  changes  significantly  in  the 
direction  of  lift  stall  delay  with  increasing  sweep,  consistent  with  Fig.  3 
of  Ref.  5  which  is  based  on  the  data  of  Ref.  6.)  A  comparison  of  the  steady 
state  stall  angle  with  the  zero  frequency  intercept  of  the  least  squares 
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vortex  inception  line  shows  a  two  degree  difference  at  Mc-  =  0.3,  and  a 
reasonably  good  correlation  at  Mc  =  0.4.  These  zero  frequency  and  steady 
state  values  are  summarized  in  Table  3. 


A  careful  study  of  Figs.  22  and  23  reveals  little  or  no  effect  of  mean 
angle  of  attack  on  the  incipient  vortex  angle.  Specifically,  at  any  constant 
value  of  kc  there  is  no  discernable  trend  of  cty  with  c^;  there  is  only 
scatter.  This  is  an  observation  that  is  directly  concerned  with  the  effect 
of  pitch  rate  on  Oy,  which  is  addressed  next. 


A  case  to  be  considered  in  detail  is  represented  by  the  solid  symbols  at 
kc  =  0.125  in  the  lower  left  panel  of  Fig.  22  for  Mc  =  0.3,  a  =  10  deg,  and 
A  *  30  deg.  In  this  instance  the  vortex  inception  angle  is  very  nearly  16 
deg  for  all  three  values  of  a^.  This  is  demonstrated  in  Fig.  24  which 
depicts  three  motion  time  histories  for  =  9,  12,  and  15  deg  versus 
dimensionless  time,  t.  A  horizontal  line  at  ay  =  16  deg  represents  the 
observed  insensitivity  of  the  vortex  inception  angle  to  changes  in  c^,  and 
leads  directly  to  a  consideration  of  pitch  rate  effects.  For  a  sinusoidal 
motion  involving  instantaneous  angle  of  attack,  a,  mean  angle  of  attack, 

0^,  and  amplitude  of  motion,  a  as  given  by  Eq .  (6),  the  instantaneous 
dimensional  pitch  rate  is  given  by 


a 


ua  cos  cot  =  oja 


-sin  ut 


(12) 


A  combination  of  Eqs.  (6)  and  (12)  yields 

a=  u>  aJol  ^  -  (a-a^)^  (13) 

which  can  be  converted  to  the  dimensionless  pitch  rate.  A,  by  the  formula 
(cf .  Ref.  7) 

A  °  2a~Jfc  “  kc  Va  2-  U-c^)^-  (14) 

The  values  of  A  for  the  intersections  of  the  ay  *  16  deg  line  with  the 
three  motion  curves  in  Fig.  24  are  noted  on  the  figure  and  are  seen  to  vary 
from  0.0156  to  0.0216.  (The  minor  changes  in  the  actual  measured  values  of 
ay  at  each  kc  in  Figs.  22  and  23  will  result  in  some  scatter  but  will  not 
materially  affect  the  conclusions  reached.  ) 

A  direct  examination  of  the  effect  of  mean  angle  of  attack  variations  on 
ay  with  reduced  frequency  effects  excluded  is  afforded  by  a  return  to 
Eq.  (11).  Recalling  that  this  equation  is  a  linear  least  squares  fit,  the 
scatter  of  individual  data  points  may  be  quantitatively  evaluated  by 
rewriting  the  equation  in  the  form 

e  -  ciy-  aVQ-  mkc  (15) 
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where  e  =  0  for  an  exact  linear  fit.  Values  of  kc,  tq>  ay  (from  Eq.  (10)), 
e  (from  Eq.  (15)),  and  Ay  (from  Eq.  (14))  with  a  _  ay)  are  listed  in 
Table  4,  and  e  is  plotted  versus  Ay  in  Fig.  25  for  all  vortex  inception 
events.  This  figure  emphasizes  the  irregularity  of  the  scatter,  and  leads 
directly  to  the  observation  that  the  vortex  inception  angle  is  independent  of 
the  pitch  rate  when  kc  is  fixed.  With  the  exception  of  a  few  isolated 
points  the  error,  e,  is  generally  bounded  by  |e |  <  0.5  deg,  and  never  exceeds 
i  1.0  deg. 

A  question  must  be  raised  at  this  point  on  the  possible  relationship  of 
the  vortex  inception  angle,  ay,  to  the  various  dynamic  stall  angle  concepts 
that  have  been  cited  in  the  literature.  As  far  back  as  Halfman's  1951  report 
(Ref.  8)  and  up  to  the  present  time,  several  authors  have  shown  that  the 
dynamic  stall  angle  increases  with  both  increasing  reduced  frequency  and 
pitch  rate.  One  of  the  more  complete  examinations  is  that  of  Gorraont  (Ref. 
9),  based  on  data  and  analyses  from  previous  work  (Refs.  10,  11).  Here  it  is 
shown  that  the  dynamic  moment  stall  angle,  based  on  3 C^(dynaraic)/3a  =  0, 
varies  linearly  with  the  square  root  of  the  pitch  rate,  and,  for  the  thick 
V23010-1.58  airfoil  examined,  approaches  the  steady  state  stall  angle  as  the 
zero  pitch  rate  limit  is  reached. 

To  permit  a  comparison  of  Gormont ' s  linear  fit  of  the  Ref.  10  dynamic 
stall  angle  results  with  the  current  data,  the  steady-state  stall  angle  was 
removed  from  the  vortex  inception  angle  results  of  the  present  experiment 
(av-  agg,  from  Tables  3  and  4),  and  then  plotted  versus  /X  in  Fig.  26.  The 
Mc  =  0.4  linear  fit  to  Gormont's  thick  airfoil  moment  stall  results  (lower 
left  panel  of  Fig.  11  in  Ref.  9)  are  shown  plotted  as  dashed  lines  in  the  two 
Mc  =  0.4  panels  of  Fig.  26.  Although  the*e  is  considerable  scatter  of  the 
data,  a  linear  trend  is  discernable  for  a  =  10  deg  (lower  two  panels). 
Gormont's  linear  fit  for  moment  stall  angle  starts  from  the  origin  and  lies 
roughly  parallel  to,  and  significantly  above,  the  ay  data  in  the  lower 
right  panel.  If  these  data  were  linearly  extrapolated  to  zero  pitch  rate, 
the  intercept  value  would  lie  approximately  1.5  deg  below  the  steady  state 
stall  angle,  which  is  in  direct  contradiction  to  the  result  shown  in  Table  3 
and  in  the  lower  right  panel  of  Fig.  22.  This  can  be  resolved  by  noting  that 
a  linear  behavior  of  the  data  is  observed  for  ay-  agg  plotted  versus  /X  in 
Fig.  26  or  versus  A  in  Fig.  27.  In  the  latter,  only  the  cases  for  a  *  10  deg 
are  repeated.  The  data  are  now  more  in  keeping  with  the  results  in  Fig.  22, 
and  the  original  Gormont  linear  data  fit  is  now  a  parabola.  The  apparent 
inconsistency  lies  in  the  lack  of  data  at  low  pitch  rates  (or  reduced 
frequency),  and  the  only  clear  result  from  the  Mc  "  0.4  evidence  is  that 
leading  edge  vortex  inception  precedes  moment  stall,  assuming  that  both  the 
V23010-1.58  airfoil  data  of  Ref.  9  and  the  current  NACA  0012  airfoil  data  can 
be  validly  compared  in  this  manner. 
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DISCUSSION  OF  NORMALIZATION  PROCEDURE 


Generally,  the  purpose  of  normalization  is  to  consolidate  data  to  a 
common  reference.  This  is  done  to  highlight  the  nonsimilar  aspects  of  the 
data  from  case  to  case  and  thus  to  permit  a  more  effective  analysis  of  their 
dependence  on  the  various  test  parameters.  In  well  behaved  small  disturbance 
problems,  for  example,  the  normalization  process  coalesces  the  data  to  a 
common  curve  so  that  they  become  amenable  to  a  general  qualitative  analysis. 
In  the  present  investigation,  however,  the  fluctuations  in  the  measured  pres¬ 
sures  exceed  what  is  normally  considered  "small  disturbance"  levels.  In 
realization  of  this  and  the  fact  that  much  of  the  data  was  obtained  within 
the  stalled  flow  regime,  the  following  discussion  attempts  to  identify  the 
limitations  of  the  normalization  procedure  used  in  this  report. 

An  example  of  how  well  the  normalization  procedure  works  below  stall 
with  the  current  data  base  is  depicted  in  Fig.  28a  in  which  unsteady  data 
obtained  from  the  30  deg  swept  wing  at  a  =  8  deg,  =  0  deg,  Mc  =  0.40,  and 
f  =  8  Hz  (kc  =  ncf/Mca0O  =  0.076)  are  compared  with  data  similarly 
obtained  at  Mc  =  0.30,  and  f  =  6  Hz  (also  at  kc  =  0.076).  This  example 
shows  that  the  reduced  frequency  normalization  procedure  successfully 
coalesces  the  linear  part  (first  harmonic  component)  of  the  chordwise 
response  distribution  to  a  common  curve  at  low  aerodynamic  load  (aM  =  0 
deg).  Conversely,  Fig.  28b  shows  some  deterioration  in  similarity  when 
results  for  the  same  combinations  of  Mach  number  and  frequency  are  compared 
at  aM  =  9  deg.  This  result  indicates  that  reduced  frequency  similarity 
deteriorates  at  high  angles  of  attack. 

The  data  set  reported  on  herein  and  in  Ref.  1  includes  effects  due  to 
wing  sweep  angle.  Sweep  was  accounted  for  in  the  normalization  process  by 
utilizing  the  component  of  freestream  velocity  normal  to  the  span  as  the 
normalizing  velocity.  This  simple  approach  is  the  same  that  has  been  tra¬ 
ditionally  used  in  the  analysis  of  steady-state  data.  An  example  of  the 
effectiveness  of  this  approach  to  swept  wing  normalization  in  similarity 
analysis  is  shown  in  Fig.  29a  for  the  zero  mean  angle  of  attack  case  at 
f  =  8  Hz,  a  =  8  deg,  and  Mc  =  0.40.  Figure  29b  shows  the  same  comparison  but 
at  <*M  =  9  deg.  As  in  Fig.  28b,  it  is  seen  that  the  higher  mean  angle  of 
attack  causes  a  deteriotat ion  in  the  sweep  similarity  between  normalized 
data. 


On  the  basis  of  the  examples  shown  in  Figs.  28  and  29,  it  would  appear 
that  the  present  normalization  procedure  is  generally  effective  in  high¬ 
lighting  the  nonsimilar  aspects  of  the  data,  and  is  thus  useful  in  deter¬ 
mining  the  effect  of  the  various  test  parameters  on  the  oscillatory  stall 
phenomenon . 
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It  is  useful  to  further  examine  this  normalization  scheme  by  using  it  to 
make  a  preliminary  evaluation  of  the  effects  of  f,  Mc ,  and  A  on  the  com¬ 
plete  unsteady  chordwise  response  of  the  upper  surface  at  zero  mean  angle  of 
attack.  Figures  30  through  32,  respectively,  show  the  effect  of  f,  M  ,  and 
A  at  the  lowest,  mean,  and  highest  incidence  points  of  the  motion  cycle.  For 
the  range  of  parameters  tested,  these  figures  show  that  the  general  shape  and 
magnitude  of  the  unsteady  chordwise  pressure  distribution  of  the  zero  mean 
angle  of  attack  configuration  are  in  good  agreement  at  the  moderate  pitching 
amplitude  of  a  =  8  deg. 

The  effectiveness  of  the  current  normalization  procedure  in  coalescing 
the  pressure  coefficient  distributions  helps  emphasize  the  differences  in  the 
leading  edge  region  results.  These  differences  are  especially  notable  in 
Fig.  31  where  the  effect  of  Mach  number  is  shown.  Since  conclusions  about 
the  behavior  of  oscillatory  pressures  cannot  be  made  on  the  basis  of  chord- 
wise  pressure  distributions  alone,  individual  pressure  time  histories  must 
also  be  examined  to  obtain  information  about  the  amplitude  and  phase  behavior 
over  the  entire  cycle.  An  illustration  of  this  is  shown  in  Fig.  33  where 
each  column  contains  cyclically  averaged  time  history  traces  of  the  four 
leading  edge  region  pressure  coefficients  from  the  same  cases  that  were  used 
to  produce  the  chordwise  plot  comparisons  of  Figs.  30  through  32.  From  left 
to  right  columns,  the  results  in  Fig.  33  show  the  effects  of  f  (corresponding 
to  Fig.  30),  Mc  (corresponding  to  Fig.  31)  and  A  (corresponding  to  Fig. 

32). 


Figure  33  shows  that  the  effect  of  increasing  f  on  the  pressure  coeffi¬ 
cient  time  history  is  to  increase  the  lag  of  the  response,  while  the  effect 
of  increasing  Mc  is  to  increase  the  magnitude  of  the  response  at  the  high 
and  low  pitching  angles  of  the  motion  cycle.  The  rightmost  column  of  Fig.  33 
shows  that  the  effect  of  A  on  the  pressure  coefficient  response  is  minimal 
with  the  swept  configuration  slightly  leading  the  unswept  case. 

Finally,  it  is  noted  that  the  parameter  kc  can  be  varied  either  by 
changing  f  or  Mc.  A  comparison  of  the  left  and  middle  columns  of  Fig.  33 
shows  that  the  effect  of  variation  in  kc  in  these  figures  is  dependent  on 
whether  f  or  Mc  is  varied.  Because  of  this  lack  of  uniqueness  in  the 
sensitivity  of  the  response  to  variations  in  kc ,  the  parameter  analysis  in 
this  report  will  focus  on  the  effects  of  f,  Mc,  and  A  rather  than  kc , 

Mc ,  and  A  . 
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UNSTEADY  CHORDWISE  PRESSURE  DISTRIBUTION  ANALYSIS 


Overview 

The  chordwise  pressure  distribution  analysis  in  this  section  is  limited 
to  8  degree  pitching  amplitude  data  with  special  attention  to  the  leading 
edge  region.  Three  examples  of  the  behavior  of  the  chordwise  distribution 
are  shown  in  Figs.  34  through  36  for  mean  angles  of  attack  of  0,  9,  and  15 
deg,  respectively.  Each  figure  shows  the  pressure  distributions  on  both 
surfaces  at  specific  points  of  the  motion  cycle  for  both  the  unswept  (left 
column)  and  swept  (right  column)  wing  configurations.  The  airfoil  pitch 
position  (relative  to  mean  angle  of  attack)  corresponding  to  these  distribu¬ 
tions  is  indicated  by  the  dot  on  the  adjacent  waveform  between  the  columns. 
Other  test  parameters,  common  to  these  examples,  are  Mc  =  0.30,  f  =  8  Hz 
(kc  =  0.102),  and  a  -  8  deg. 

Before  proceeding  with  the  description  of  the  results,  it  is  noted  that 
steady  state  stall  occurs  at  approximately  a  =  13  deg  for  Mc  =  0.30  and  at 
approximately  a  =  10.5  deg  for  Mc  =  0.40  (Table  3).  Therefore,  in  terms  of 
the  present  test  configurations,  dynamic  stall  penetration  occurs  at  angles 
of  attack  greater  than  the  mean  for  =  9  deg,  either  slightly  above  or 
below  the  mean  for  =  12  deg  (depending  on  Mc )  and  below  the  mean  for 
cfy  =  15  deg.  This  aerodynamic  behavior  relative  to  should  be  borne  in 
mind  during  the  discussion  of  the  examples  in  this  section. 

The  cases  shown  in  Figs.  34  through  36  are  singled  out  as  typical  of  th 
entire  data  base  for  the  range  of  mean  angles  of  attack  tested.  In  addition 
to  being  a  representative  zero  degree  mean  angle  of  attack  case,  Fig.  34  can 
also  be  used  as  a  benchmark  for  evaluating  the  quality  and  consistency  of  th 
data.  For  example,  the  upper  and  lower  surface  chordwise  pressure  distribu¬ 
tions  in  Fig.  34  appear  nearly  symmetric  with  the  motion  within  each  column. 
This  is  consistent  with  the  profile  shape  of  the  NACA  0012  airfoil.  In  par¬ 
ticular,  note  how  the  main  features  of  the  upper  and  lower  surface  responses 
are  reversed  at  mean  angle  of  attack  (a  =  0  deg)  when  the  upstroke  and  down- 
stroke  distributions  are  compared.  Similarly  note  the  reversal  of  the 
surface  distributions  at  the  maximum  and  minimum  points  of  the  motion  cycle. 
Another  test  of  consistency  is  to  compare  the  unswept  and  swept  results.  As 
expected  by  the  discussion  of  the  previous  section,  the  data  appear  nearly 
the  same  throughout  the  motion  cycle.  (It  is  further  noted  that  the  unswept 
and  swept  wing  results  were  generated  in  two  separate  tunnel  entries.) 

In  contrast,  the  stalled  flow  penetration  case  shown  in  Fig.  35  (a^  = 

9  deg)  shows  no  evidence  of  symmetry  in  the  response  about  the  mean  angle  of 
attack.  This  example  also  shows  that  sweep  similarity  is  s igni f icant ly 


altered  in  the  higher  incidence  angle  range  as  the  airfoil  oscillates  in  and 
out  of  stall.  In  this  comparison,  the  decrease  in  similarity  is  especially 
noticeable  at  peak  incidence  where  the  unswept  wing  has  experienced  a  more 
severe  collapse  of  the  suction  response  than  the  swept  configuration. 

Figure  36  (a^  =  15  deg)  shows  chordwise  distributions  that  were 
generated  mostly  within  the  stalled  flow  region.  Here  it  is  seen  that  the 
differences  between  the  unswept  and  swept  wing  results  are  not  as  severe, 
giving  the  appearance  that  some  sweep  similarity  is  regained  when  the  airfoil 
is  in  deep  stall.  However,  as  will  be  shown  in  later  paragraphs  of  this 
section,  this  is  not  the  case  since  there  remains  substantial  disagreement, 
especially  in  the  forwardmost  2  percent  region  of  the  chord,  during  the 
entire  motion  cycle.  On  the  basis  of  Figs.  35  and  36  it  appears  that  the 
dissimilarities  due  to  sweep  are  strongest  when  the  airfoil  load  is  in 
transition  between  an  unstalled  and  a  stalled  response. 


Effect  of  on  Normalization 

The  remaining  discussion  of  the  chordwise  pressure  data  is  focused  on 
the  sensitivity  of  the  suction  (upper)  surface  distribution  to  the  various 
test  parameters  with  occasional  emphasis  on  the  effects  of  sweep.  The  cases 
and  their  behavior  highlighted  in  this  section  are  typical  of  the  entire  data 
set . 


In  the  previous  section  it  was  shown  that  the  normalization  procedure 
used  in  the  data  reduction  successfully  coalesces  the  first  harmonic  compon¬ 
ent  of  the  zero  mean  incidence  angle  cases  to  a  common  curve  (cf.  Figs.  28 
and  29).  This  coalescence  was  further  demonstrated  in  Figs.  30  through  32 
for  the  complete  unsteady  component  of  the  total  response  at  zero  mean  angle 
of  attack.  In  these  latter  three  figures  comparisons  were  made  showing  the 
respective  effects  of  pitching  frequency,  Mach  number,  and  sweep. 

The  impact  of  mean  angle  of  attack  on  this  generalized  behavior  is  now 
examined  by  comparing  the  swept  wing  results  in  Figs.  37  (aM  =  0  deg)  and 
38  (a^  =  9  deg).  Each  figure,  in  turn,  shows  comparisons  between  two  test 
cases  in  which  both  the  frequency  and  Mach  number  are  varied  simultaneously 
while  maintaining  a  constant  reduced  frequency  of  kc  *  0.076.  The  varia¬ 
tions  in  frequency  and  Mach  number  between  cases  are  the  same  in  each 
figure. 

There  is  nearly  perfect  agreement  between  the  two  cases  over  the  entire 
cycle  at  aM  *  0  deg  (Fig.  37)  confirming  the  validity  of  reduced  frequency 
scaling  at  low  load.  However,  when  a  similar  comparison  is  made  at  aM  = 

9  deg  (Fig.  38),  the  two  cases  are  in  disagreement.  This  deterioration  in 
generalized  behavior  occurs  over  the  entire  high  angle  of  attack  range  of  the 
motion  cycle  with  maximum  disagreement  at  peak  a  where  the  unsteady  separ¬ 
ation  effects  spread  over  the  entire  chord.  Conversely,  it  is  noted  that 
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similarity  in  the  chordwise  response  is  regained  at  the  minimum  point  of  the 
cycle  where  a  =  1.0  deg  (Fig.  38).  The  disagreement  in  the  chordwise 
responses  over  most  of  the  cycle  is  due  to  the  difference  in  Mach  number 
between  the  two  cases  as  illustrated  in  the  following  discussion. 

The  sensitivity  of  the  data  at  =  9  deg  to  variations  in  the  reduced 

frequency  is  examined  in  the  next  two  figures  for  both  the  unswept  (left 

column)  and  swept  (right  column)  wings  at  the  mean,  maximum,  and  minimum 

points  of  the  motion  cycle.  The  objective  here  is  to  show  the  influence  of 

kc  on  the  data  by  separately  varying  the  pitching  frequency  (Fig.  39)  and 

Mach  number  (Fig.  40).  In  Fig.  39,  the  variation  in  kc  (due  to  f)  is  three 

times  that  of  Fig.  40  (due  to  Mc)  and  yet  the  degree  and  trend  of  the 

change  in  the  chordwise  response  between  cases  in  each  figure  is  about  the 

same  for  both  the  unswept  and  swept  configurations.  These  figures  show  that 

the  trend  of  the  pressure  response  distribution,  although  different  in 

magnitude,  is  generally  consistent  with  the  variation  in  kc  (due  to  changes 

in  either  f  or  M  )  and  verify  the  dominant  influence  of  M  on  k  . 

c  J  c  c 

A  further  inspection  of  Figs.  39  and  40  also  shows  that  although  the 
influence  of  kc  (due  to  changes  in  either  f  or  M  )  is  generally  the  same 
between  unswept  and  swept  results,  there  are  regions  of  the  motion  cycle 
where  kc  has  the  opposite  effect.  (For  example,  near  peak  angle  of  attack 
in  Fig.  40,  an  increase  in  kc  corresponds  to  a  decrease  in  suction  for  the 
unswept  response  and  an  increase  in  suction  for  the  swept  response  in  the 
forward  10  percent  region  of  the  chord.) 

This  difference  in  the  effect  of  kc  on  the  response  at  spec':-  -  poin*? 
of  the  cycle  is  partly  due  to  the  difference  in  the  phasing  of  W  .snsteaJy 
load  response  relative  to  the  pitching  motion  between  the  unswept  and  swept 
conf igurat ions .  As  a  result,  the  effect  of  kc  shows  up  differently  at 
various  points  of  the  cycle,  especially  in  the  transition  region  (between 
stalled  and  unstalled  flow)  where  the  airfoil  surface  pressure  response 
appears  to  be  most  sensitive  to  parameter  variations.  Examples  illustrating 
the  sensitivity  of  the  phasing  of  the  leading  edge  region  pressure  response 
are  presented  next. 


Leading  Edge  Region  Time  Histories 

An  important  aspect  in  the  parameteric  evaluation  of  chordwise  load 
behavior  is  the  analysis  of  the  pressure  time  histories  at  various  points 
along  the  chord.  The  objective  here  is  to  show  representative  examples  from 
the  data  base  and  to  highlight  some  of  the  main  features  in  the  behavior  of 
the  pressure  time  histories  to  variations  in  the  test  parameters.  The  time 
histories  shown  in  Figs.  41  through  43  each  contain  three-harmonic  component 
reconstructions  of  the  leading  edge  region  responses  at  the  0.4,  1.0,  and  1.9 
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percent  chord  positions  of  the  suction  surface.  (The  choice  of  three-har¬ 
monic  component  plots  rather  than  the  total  response  is  to  keep  the  para¬ 
metric  analysis  relatively  simple  while  not  giving  up  any  essential 
features.)  Each  figure  is  arranged  in  3  columns  corresponding  to  the  0,  9, 
and  15  degree  mean  incidence  angle  configurations.  Each  column  is  comprised 
of  the  blade  pitching  time  history  at  the  top  followed  by  the  three  forward- 
most  suction  surface  time  histories  at  x  =  0.004,  0.010,  and  0.019. 

Figure  41  shows  the  effect  of  frequency  at  Mc  =  0.30,  A  =  0  deg,  and 
a  =  8  deg.  In  this  figure,  it  is  seen  that  a  frequency  increase  from  4 
to  8  Hz  (or  kc  increasing  from  0.051  to  0.102)  has  only  a  minor  impact, 
showing  up  mainly  as  a  slight  increase  in  the  phase  lag  (right  shift)  of  the 

responses  at  =  0  deg  and  9  deg.  At  =  1 5  deg,  however,  this 

increase  in  frequency  causes  a  substantial  increase  in  the  phase  lag  of  the 
response  during  the  upstroke  as  well  as  a  significant  dimunition  in  the 
response  level  during  the  downstroke.  Figure  44  is  a  companion  to  Fig.  41 
where  the  9  deg  (left  column)  and  15  deg  (right  column)  mean  incidence  angle 
time  histories  from  the  leading  edge  region  are  repeated  with  additional  time 
histories  further  down  along  the  chord  at  x  =  0.045,  0.073,  0.114,  and  0.268. 
This  figure  shows  little  variation  in  the  relative  phasing  throughout  the 

forwardmost  quarterchord  region.  This  is  most  easily  seen  during  the  first 

quarter  cycle  of  the  =  15  deg  configuration. 

Figures  42  and  45  show  the  effect  of  Mach  number  at  f  =  8  Hz,  =  0  deg, 
and  a  =  8  deg  using  the  same  comparison  approach  as  in  Figs.  41  and  44.  Here 
Fig.  42  shows  that  an  increase  in  Mach  number  has  a  greater  effect  on  the 
leading  edge  region  response  than  does  the  variation  in  the  frequency  (cf. 
Fig.  41).  This  sensitivity  of  the  leading  edge  region  response  to  the  Mach 
number  is  especially  significant  at  aM  =  9  deg  and  to  a  lesser  extent  at 
aM  =  15  deg.  In  particular,  it  is  noted  that  the  increase  in  Mach  number 
lowers  the  stall  penetration  angle  and  reduces  the  sharpness  of  the  dynamic 
stall  event.  Figure  45  shows  that  this  sensitivity  of  the  aerodynamic 
response  to  M£  persists  at  least  to  x  5  0.268. 

Examples  of  the  effect  of  sweep  on  the  time  history  behavior  of  the 
chordwise  response  are  shown  next  in  Figs.  43  and  46  at  f  =  8  Hz ,  Mc  = 

0.30,  and  a  =  8  deg.  The  zero  degree  mean  angle  of  attack  results  in  Fig.  43 
show  that  sweep  effects  are  essentially  negligible  throughout  the  cycle  all 
way  to  the  forwardmost  station  at  x  “  0.004.  Conversely,  sweep  signficantly 
alters  the  leading  edge  region  response  when  =  9  and  15  deg.  In  the 
9  deg  mean  angle  of  attack  case,  sweep  effects  remain  fairly  significant 
through  the  quarterchord  region  while  the  15  deg  mean  angle  of  attack  case 
shows  a  substantial  decline  in  the  influence  of  sweep  beyond  x  “  0.045.  The 
most  notable  effect  of  sweep  on  these  results  is  in  the  increased  phase  lag 
of  the  load  and  the  associated  delay  in  the  stall  event.  This  is  especially 
clear  when  aM  "  9  deg.  In  particular,  note  how  the  peak  response  during 
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the  upstroke  is  delayed  to  later  in  the  cycle  for  A  =  30  deg  (dashed  line) 
relative  to  A  =  0  deg  (solid  line). 


Parameter  Analysis  of  Chordwise  Pressures 

The  effect  of  mean  angle  of  attack  on  the  suction  surface  pressure 
response  and  on  the  influence  of  sweep  at  specific  angles  of  attack  is 
examined  next.  As  examples,  Figs.  47  and  48  show  this  effect  at  a  =  7.4  deg 
and  16.6  deg  in  the  upstroke  portion  of  the  cycle.  Each  figure  shows  results 
at  Mc  =  0.30  (left  column)  and  at  Mc  =  0.40  (right  column)  along  with 
accompanying  simple  waveforms  indicating  the  point  in  the  motion  cycle  where 
the  chordwise  distributions  were  taken.  Figure  47  contains  results  taken 
from  the  portion  of  the  cycle  where  the  unsteady  aerodynamic  activity  is  not 
in  transition  and  is  within  the  potential  flow  range  of  the  motion.  It  is 
seen  that  the  unswept  and  swept  wing  chordwise  distributions  in  this  part  of 
the  cycle  are  each  nearly  similar  over  the  range  of  mean  angle  of  attack 
shown  at  Mach  numbers  of  0.30  and  0.40.  This  behavior,  although  not 
completely  generalized,  is  encouraging  from  a  data  utilization  standpoint  in 
that  this  response  distribution  appears  applicable  (at  this  point  of  the 
motion  cycle)  over  the  entire  range  of  mean  incidence  angles  between  = 

9  deg  and  =  15  deg.  With  the  exception  of  the  upper  left  distributions, 
Fig.  48  shows  that  this  similarity  of  the  unswept  and  swept  wing  responses 
with  respect  to  mean  angle  of  attack  tends  to  persist  into  the  separated  flow 
region.  In  particular,  the  chordwise  distributions  for  each  Mach  number  in 
Fig.  48  are  very  similar  at  =  12  deg  and  15  deg  while  noticeably  differ¬ 
ent  at  *  9  deg  where  the  angle  of  attack  value  used  in  this  comparison 
is  near  its  peak.  (It  is  again  pointed  out  here  that  the  jaggedness  in  the 
leading  edge  region  response  is  an  anomaly  of  the  data  which  was  explained 
earlier.)  It  is  also  noted  in  this  figure  that  the  effects  of  Mach  number 
and  sweep  are  substantial  at  this  high  value  of  angle  of  attack  (16.6  deg), 
while  the  effect  of  Mach  number  at  a  =  7.4  deg  (cf.  Fig.  47)  is  nearly 
impercept ible . 

A  further  analysis  of  the  Mc  =  0.30  results  corresponding  to  the 
previous  two  figures  is  presented  in  Figs.  49  and  50.  As  in  Figs.  47  and  48, 
these  figures  show  the  effects  of  sweep  and  mean  angle  of  attack  on  the 
chordwise  distribution  at  various  points  of  the  motion  cycle.  Figure  49 
shows  the  chordwise  distributions  at  instantaneous  angles  of  attack  of  a  =  9, 
11,  13,  and  15  deg  during  upstroke  (in  columns  from  left  to  right)  thus  fill¬ 
ing  the  ga  between  the  results  of  Figs.  47  and  48  at  7.4  deg  and  16.6  deg 
incidence,  respectively.  It  is  seen  that  the  effect  of  varying  on  the 
shape  of  the  distribution  (within  each  column)  is  relatively  small  over  this 
range  of  angle  of  attack  (within  each  row)  for  both  the  unswept  and  swept 
cases,  implying  that  significant  stall  penetration  has  not  yet  taken  place. 

It  is  also  noted  that  the  effect  of  sweep  in  reducing  suction  is  consistent 
throughout  this  portion  of  the  motion  cycle. 
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A  similar  analysis  is  presented  in  Fig.  50  for  the  downstroke  portion  of 
the  cycle  over  the  same  angle  of  attack  range  (where  the  plots  are  arranged 
following  the  motion  from  left  to  right).  In  this  region  the  effects  of 
sweep  and  mean  angle  of  attack  are  each  more  significant  than  in  Fig.  49. 

This  is  especially  evident  at  a  =  15  deg  (extreme  left  column)  where  flow 
reattachment  is  still  underway.  The  sensitivity  of  the  chordwise  response  to 
changes  in  mean  angle  of  attack  during  the  downstroke  portion  of  the  motion 
cycle  occurs  mainly  because  the  aerodynamic  activity  is  in  varying  stages  of 
recovery  from  separation.  This  variation  in  the  level  of  aerodynamic  activ¬ 
ity  depends  both  on  how  deeply  the  stalled  flow  region  was  penetrated  by  the 
airfoil  and  on  where  the  airfoil  is  in  the  motion  cycle  relative  to  peak 
incidence. 


Effect  of  Sweep  on  Peak  Angle  of  Attack  Response 

The  effect  of  sweep  on  the  peak  angle  of  attack  region  response  is  now 
examined  for  aM  =  9,  12  and  15  deg.  Figures  51  and  52  show  the  sensitivity 
of  the  chordwise  response  to  sweep  in  0.5  degree  increments  within  1.5  deg  of 
peak  angle  for  each  mean  angle  of  attack  case.  Figure  51  shows  the  upstroke 
portion  of  the  study  while  Fig.  52  shows  the  downstroke  portion. 

The  effect  of  sweep  on  the  9  degree  mean  angle  of  attack  configuration 
is  examined  first  (see  top  row  of  Figs.  51  and  52).  At  a  =  15.5  deg  in  the 
top  row  in  Fig.  51,  the  effect  of  sweep  is  still  moderate  as  the  airfoil 
motion  enters  the  stall  transition  region  of  the  cycle.  A  half  degree  later 
at  a  *  16.0  deg,  there  is  beginning  evidence  of  collapse  in  suction  for  the 
unswept  wing.  This  collapse  progresses  rapidly  to  maximum  severity  at  peak 
angle  while  the  loss  in  suction  for  the  swept  wing  occurs  in  a  more  orderly 
fashion  and  does  not  reach  the  same  intensity.  As  a  result  of  this  differ¬ 
ence  in  the  rate  and  magnitude  of  suction  collapse,  the  unswept  and  swept 
wing  chordwise  load  distributions  are  signficantly  different  at  peak  angle  of 
attack.  Going  on  to  the  top  row  in  Fig.  52  where  the  downstroke  portion  is 
presented,  the  change  in  both  the  unswept  and  swept  wing  responses  becomes 
leas  severe  with  the  change  in  angle  of  attack  and  the  chordwise  distribu¬ 
tions  appear  to  once  again  coalesce  as  a  reaches  15.5  degrees. 

The  results  for  the  12  degree  mean  angle  of  attack  configuration  are 
shown  in  the  middle  row  of  both  figures.  In  this  case  the  motion  is  well 
into  the  stall  region  throughout  the  range  of  angles  shown  and  the  pressure 
distribution  is  still  exhibiting  the  after-effects  of  transition  between 
unstalled  and  stalled  flow.  This  behavior  is  made  more  apparent  by  the 
extensive  disagreement  between  the  unswept  and  swept  airfoil  pressure  distri¬ 
butions  throughout  the  upstroke  portion  of  the  motion  near  peak  angle  of 
attack  (Fig.  51).  This  sweep-induced  difference  becomes  less  severe  as  the 
motion  passes  through  the  peak  and  into  the  downstroke  region  (Fig,  52). 
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However,  the  strong  disagreement  at  the  leading  edge  persists  throughout  this 
range  of  angles.  As  shown  earlier  in  the  discussion  of  response  time 
histories,  this  behavior  is  partly  a  consequence  of  the  difference  in  the 
phasing  of  the  recovery  process  between  the  unswept  and  swept  wing  responses 
relative  to  the  airfoil  pitching  motion. 

The  behavior  of  the  15  degree  mean  angle  of  attack  configuration  is 
shown  in  the  bottom  row  of  Figs.  51  and  52.  In  this  case,  the  chordwise 
distributions  all  occur  very  deep  in  the  stall  region,  where  flow  separation 
is  fully  achieved,  and  are  thus  relatively  stable  to  small  changes  in  angle. 
At  a  =  21.5  deg  in  Fig.  51,  it  is  seen  that  the  transition  into  fully 
separated  flow  is  nearly  completed.  Once  the  transition  is  fully  achieved, 
the  unswept  and  swept  distributions  appear  to  coalesce  and  as  the  motion 
passes  through  peak  angle  of  attack,  the  load  distributions  appear  less 
dependent  on  sweep.  It  should  be  noted,  however,  that  this  latter 
observation  can  be  misleading  if  it  is  not  combined  with  a  study  of  the 
corresponding  time  histories  of  Fig.  43  where  it  is  seen  that  the  close 
agreement  between  unswept  and  swept  results  is  only  transitory  during  the 
downstroke  at  T  =  0.33  (right  column). 


DISCUSSION  OF  INTEGRATED  LOAD  RESULTS 


The  integrated  load  results  (lift,  drag,  and  moment)  reported  in  Ref.  1 
were  computed  from  the  chordwise  response  data  of  this  report.  Hysteresis 
loop  comparisons  were  made  therein  to  show  the  effect  of  sweep  on  the 
unsteady  lift,  drag,  and  moment  components  of  the  total  force  response.  A 
particular  feature  that  remained  unexplained  in  the  results  is  the  substan¬ 
tial  displacement  in  the  lift  curve  below  stall  between  the  unswept  and  swept 
wing  results,  especially  at  high  angles  of  attack.  A  perspective  on  this 
displacement  can  be  obtained  from  Figs.  53  and  54.  Figure  53  shows  that  the 
normalized  steady-state  lift  curves  of  the  unswept  and  swept  NACA  0012  air¬ 
foils  are  nearly  coalesced  below  stall.  This  is  not  the  case  when  viewing 
the  unsteady  comparisons  in  Fig.  54.  Here,  it  is  seen  that  the  upstroke 
portions  of  the  unswept  lift  responses  below  stall  are  significantly  and 
increasingly  displaced  from  one  another  as  increases.  (Figure  54  is  a 
reproduction  of  results  from  Fig.  34  of  Ref.  1  showing  the  effect  of  sweep 
on  the  lift  response  for  four  mean  angles  of  attack  tested  at  Mc  =  0.40, 
f  =  Hz ,  and  a  =  8  deg.) 

It  has  been  found  that  this  displacement  of  the  lift  reponse  is  mainly 
associated  with  the  first  harmonic  component.  This  was  discovered  serendipi- 
tously  while  attempting  to  superimpose  the  unsteady  parts  of  the  unswept  and 
swept  wing  lift  responses.  In  order  to  achieve  this  superposition  effec¬ 
tively,  the  time  mean  components  of  the  lift  responses  were  removed  in  an 
attempt  to  bring  the  loops  to  closer  alignment  than  shown  in  Fig.  54.  The 
result  of  doing  this  is  shown  in  the  top  row  of  Fig.  55  for  the  15  degree 
mean  angle  of  attack  case  where  plot  55a  is  repeated  from  Fig.  54  and  plot 
55b  is  stripped  of  the  time  mean  component.  The  resultant  lack  of  change  in 
the  relative  positioning  of  these  loops  indicates  that,  for  this  particular 
comparison,  the  entire  displacement  is  a  sweep-induced  unsteady  phenomenon. 

A  harmonic  analysis  of  these  loops  was  performed  to  determine  the  aspect  of 
the  data  that  accounts  for  this  displacement.  It  was  found  that  the  ultimate 
shape,  displacement,  and  rotation  of  the  lift  loops  of  Fig.  55a  are  primarily 
governed  by  the  first  harmonic  component  ellipse.  This  is  clearly  illus¬ 
trated  by  the  first  harmonic  component  comparison  in  plot  55c  where  the  mid¬ 
portion  of  the  AC^  upstroke  displacement  of  the  complete  response  (plot 
55a)  i 8  accounted  for.  This  first  harmonic  result  shows  that  the  primary 
influence  of  sweep  on  the  total  lift  response  is  through  the  phasing  between 
the  first  harmonic  component  of  the  lift  response  and  the  blade  motion,  and 
through  the  ratio  of  the  first  harmonic  component  amplitude  of  the  response 
to  the  airfoil  pitching  amplitude.  The  remaining  comparisons  in  Fig.  55 
(plots  d  through  h)  show  how  the  addition  of  higher  harmonic  terms  contrib¬ 
utes  in  evolving  the  complete  unsteady  response  without  significant  influence 
on  the  AC^  displacement  at  the  midpoint  of  the  upstroke.  From  this 
example,  it  is  seen  that  the  main  importance  of  the  higher  harmonic  terms  is 
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in  shaping  the  loops  rather  than  affecting  their  orientation.  A  comparison 
of  plot  55d  and  plot  55b  shows  that  most  of  the  remaining  features  of  the 
loops  are  due  to  the  second  harmonic  term. 

The  analysis  of  the  moment  responses  for  the  same  unswept  and  swept 
cases  shown  in  Fig.  55  is  presented  in  Fig.  56.  Unlike  the  lift  result,  it 
is  seen  that  the  removal  of  the  nonoscillatory  terras  from  the  total  values 
tends  to  bring  the  unswept  and  swept  responses  to  closer  proximity  but  with 
some  displacement  remaining.  As  in  the  case  of  the  lift  result  the  subse¬ 
quent  harmonic  analysis  of  the  moment  response  shows  that  the  remaining  dis¬ 
placement  is  due  to  the  influence  of  sweep  on  the  unsteady  portion  of  the 
response  and  that  the  ultimate  orientations  are  again  largely  determined  by 
the  first  harmonic  ellipses.  In  this  case  the  influence  of  the  second 
harmonic  term  (plot  56d)  on  the  final  shape  is  clearly  evident. 

Several  additional  sweep  comparisons  were  analyzed  and  variations  were 
observed  in  the  source  (unsteady  or  time  mean)  of  the  loop  displacement  from 
case  to  case.  An  example  that  illustrates  this  variation  in  the  sweep  influ¬ 
ence  (relative  to  that  shown  in  Figs.  55  and  56)  is  presented  next  in  Figs. 

57  and  58. 

The  unswept  and  swept  cases  shown  in  these  figures  are  the  same  in  every 
respect  to  Figs.  55  and  56,  except  that  =  9  deg.  The  results  of  the 
harmonic  analysis  of  the  lift  in  Fig.  57  is  basically  the  •  ame  as  before; 
however,  Fig.  58  shows  that,  in  this  case,  the  displacement  of  the  moment 
response  is  entirely  accounted  for  by  the  nonoscillatory  terms.  As  a  result, 
note  how  the  first  harmonic  ellipses  are  nearly  coincident. 

Taken  together,  the  harmonic  analysis  presented  in  Figs.  55  through  58 
show  that  the  loop  displacements  range  from  being  entirely  dependent  on  the 
unsteady  influence  (cf.  Fig.  55)  to  entirely  dependent  on  the  time  mean 
influence  (cf.  Fig.  58).  Therefore,  what  once  appeared  to  be  anomalous 
behavior  (because  of  the  inconsistent  result  obtained  by  removal  of  the  time 
mean  component  from  the  total  response)  has  been  shown  to  be  amenable  to 
analysis . 
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SUMMARY  OF  RESULTS 


An  analysis  was  performed  on  unsteady  aerodynamic  pressure  data  from  a 
NACA  0012  tunnel-spanning  wing  oscillating  in  the  unswept  and  swept  configu¬ 
ration.  The  range  of  parameters  tested  were:  sweep,  0  and  30  degrees; 
pitching  amplitude,  8  and  10  degrees;  mean  angle  of  attack,  0  to  15  degrees; 
Mach  number  normal  to  the  wing  span,  0.30  and  0.40;  pitching  frequency,  4  to 
11  Hertz.  The  corresponding  reduced  frequency  range  of  the  data  was  .038  to 
.125. 


The  main  subject  areas  of  this  report  are:  the  analysis  of  the  chord- 
wise  wave  propagation  speed,  the  analysis  of  the  upper  surface  chordwise 
pressure  distribution,  and  a  brief  harmonic  analysis  of  the  integrated  lift 
response.  The  results  of  this  study  are  summarized  below.  Items  with 
asterisks  are  regarded  as  primary  findings. 


Pressure  Wave  Analysis 

•  Mean  angle  of  attack  has  little  or  no  effect  on  wave  speed,  which  in 
all  cases  increases  uniformly  with  reduced  frequency  (Figs.  16  through 
20).  Motion  amplitude  also  has  little  or  no  effect  on  wave  speed  (Fig. 
16).  Mach  number  has  a  small  but  consistent  effect  on  wave  speed,  with 
higher  Mach  numbers  yielding  slightly  smaller  wave  speeds  (Fig.  17). 

*•  Sweep  angle  has  a  dominant  effect  on  wave  speed.  The  overall  wave 

speed  for  the  unswept  wing  is  consistently  greater  than  that  for  the 
swept  wing  by  a  factor  of  1.5  to  2.0  (Figs.  18  et  seq.).  This 
represents  a  major  failure  above  stall  for  the  cosine  law  normalization 
which  has  been  shown  to  be  consistently  valid  below  stall. 

•  Previous  observations  of  the  effect  of  increasing  sweep  to  delay  stall 
are  confirmed  in  the  present  work.  In  addition,  the  vortex  inception 
angle  increases  slowly  with  kc  (Fig.  22). 

*•  Mach  number  has  a  dominant  effect  on  vortex  inception  angle  and  static 
stall  angle,  with  a  decrease  in  both  angles  as  Mc  increases  (Fig. 

23). 

•  Vortex  inception  is  substantially  independent  of  amplitude  of  motion 
and  of  mean  angle  of  attack  (Figs.  22,  23). 

•  Local  wave  speed  differences  associated  with  sweep  are  confined  to  the 
region  of  the  blade  aft  of  the  15  percent  chord.  Wave  speeds  are 
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substantially  the  same  for  both  sweep  angles  forward  of  this  chord  loca¬ 
tion  (Figs.  19,  20).  The  variation  of  local  wave  speed  with  chordwise 
position  shows  an  initial  rapid  increase  near  the  leading  edge  followed 
by  a  leveling  off  of  wave  speed  over  the  aft  region  of  the  chord  (Fig. 
21).  Overall  wave  speed  (averaged  over  the  first  15  percent  of  the 
chord)  agrees  well  with  previous  results  from  Ref.  4  (Fig.  15). 

*•  The  vortex  inception  angle  is  independent  of  pitch  rate  when  reduced 
frequency  is  fixed  (Fig.  25). 

Chordwise  Pressures  Analysis 

•  Sweep  effects  on  the  upper  surface  pressure  distribution  are  more 
prevalent  when  the  aerodynamic  response  is  in  transition  between 
unstalled  and  stalled  conditions. 

•  The  sensitivity  of  the  pressure  responses  to  sweep  is  also  strongest 
when  the  wing  angle  of  attack  is  near  peak  value  for  mean  angles  of 
attack  of  9,  12,  and  15  degrees  (Figs.  34,  35,  and  36), 

*•  Sweep  similarity  by  means  of  the  cosine  law  (which  prevails  over  the 

angle  of  attack  range  of  i  8  degrees  about  zero  mean  angle)  deteriorates 
when  the  mean  angle  of  attack  reaches  9  degrees  and  the  motion  pene¬ 
trates  stall  (Figs.  43  and  46). 

•  The  effects  of  sweep  very  near  peak  angle  appear  stronger  on  the 
approach  to  than  on  recession  from  peak  angle  when  aM  =  9  deg  while 
the  opposite  appears  to  be  the  case  when  aM  =  12  deg.  The  effects  of 
sweep  at  aM  =  15  deg,  on  the  other  hand,  are  relatively  small  near 
peak  angle  (Figs.  51  and  52). 

*•  The  reduced  frequency  parameter  is  not  a  reliable  indicator  of  the 
unsteady  aerodynamic  condition  at  high  angles  of  attack.  This  is 
because  the  same  numerical  value  of  reduced  frequency  can  be  obtained 
for  two  different  combinations  of  Mc  and  f  in  which  the  Mc  effect  on 
the  flow  dominates  the  effect  of  f  (Figs.  37  and  38). 

Integrated  Load  Analysis 

*  •  The  main  influence  of  sweep  on  the  total  integrated  lift  and  moment 

responses  is  very  strongly  determined  by  the  first  harmonic  term.  This 
result  implies  that  the  effect  of  sweep  on  the  unsteady  response  is 
amenable  to  analysis  (Fig.  55). 
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APPENDIX 
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DISCUSSION  OF  MEASUREMENT  ERROR 


A  minor  measurement  error  was  detected  during  the  present  analysis  that 
occurs  throughout  the  data  set.  This  error  neither  affects  the  discussion 
presented  herein  nor  does  it  alter  the  conclusions  previously  reported  in 
Ref.  1  (based  on  the  same  data  set). 

The  main  evidence  of  the  error  is  the  jaggedness  that  occurs  in  both  the 
upper  surface  chordwise  pressure  distribution  plots  and  the  constant  pressure 
contours  of  the  chordwise  wave  speed  analysis.  For  the  chordwise  distribu¬ 
tion  plots,  this  anomaly  ranges  from  slight  to  relatively  large  depending  on 
the  pressure  loading  condition.  In  the  case  of  the  constant  pressure  con¬ 
tours  the  effect  is  generally  minor  and,  at  worst,  is  more  a  nuisance  than  a 
detriment  to  the  evaluation  of  the  wave  speed. 

A  close  inspection  of  the  first  harmonic  component  phase  angle  distribu¬ 
tion  along  the  chord  in  Fig.  4  shows  that  the  jaggedness  (seen  in  Fig.  5  for 
example)  is  principally  related  to  a  phasing  error  between  consecutive  pres¬ 
sure  responses  along  the  chord.  In  particular,  it  is  noted  that  the  phase 
angle  plots  near  the  leading  edge  appear  to  be  scattered  about  a  mean  trend. 
This  behavior  of  the  data  set  was  found  to  be  caused  by  a  misalignment  of  the 
tape  heads  of  the  analog  tape  recorders.  Both  the  read  and  write  functions 
were  configured  to  utilize  two  heads  during  record  and  playback  where  the  odd 
and  even  numbered  channels  (sequentially  assigned  along  the  chord)  were 
segregated  to  separate  heads. 

Under  these  circumstances,  the  time  lag  between  odd  and  even  channels 
due  to  head  misalignment  caused  an  odd-even  channel  phase  error.  This  phase 
error,  A  ,  is  related  to  the  oscillation  period,  T,  and  tape  speed,  S,  by  the 
equat ion 


=  At  =iL  =  f  it  s  f 
360°  T  S 

where  Ax  represents  the  error  displacement  between  the  heads.  A  typical 
value  of  Ay  equal  to  .003  in.  produces  a  phase  angle  error  of  5  deg  in  the 
first  harmonic  at  f  =  8  Hz. 

The  entire  data  set  was  uniformly  corrected  for  the  computations  made  in 
the  wave  speed  analysis,  but  the  corrections  were  confined  to  that  section 
only.  All  other  figures  presented  in  this  report  and  the  data  tables  in  Ref. 
3  remain  uncorrected  to  preserve  consistency  with  the  results  of  Refs.  1 
and  2 . 
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TABLE  I 


UNSTEADY  TEST  MATRIX  FOR  NACA0012  AIRFOIL 
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TABLE  2 


NOMINAL  VALUES  OF  kc  FOR  BASIC  TEST  PROGRAM 


All  values  displayed  were  run  at  least  once. 

Horizontal  arrows  on  right  column  denote  nominal  values 
used  for  unswept  runs. 

Boxed  numbers  denote  nominal  values  used  for  swept  runs 
Diagonal  arrows  indicate  matched  values  of  kc. 


TABLE  3 


Least  Squares 
intercept  at 

kc=0 

Steady  State 

Based  on 
dC^/da  =0 

Mc 

a 

A=0° 

A=30  ° 

o 

O 

II 

< 

A=30  ° 

.3 

8° 

14.4“ 

15  .0° 

12  .9° 

12  .8 6 

10° 

14.5° 

14  .5“ 

.4 

8° 

10.3° 

11  .4“ 

o 

o 

10  .4° 

10° 

9  .9° 

10.6° 
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TABLE  4 


,aM  \ 
(deg) 

kc 

T 

O 

av 

( deg) 

G 

(deg) 

Av 

9 

.0498 

.137 

15  .07 

.423 

.0045 

A=0  deg 

.0993 

.125 

14  .66 

-.202 

.0098 

.1245 

.123 

14.58 

-.391 

.0125 

M  =0.3 

12 

.0497 

.052 

14 .57 

-.077 

.0066 

.1001 

.057 

14  .80 

-.065 

.0131 

5=8  deg 

.1249 

.066 

15  .22 

.247 

.0160 

15 

.0497 

.987 

14 .35 

-.297 

.0069 

.1002 

.000 

15.00 

.134 

.0140 

.1248 

.004 

15.20 

.228 

.0174 

9 

.0499 

.094 

14.57 

-.334 

.0072 

A=0  deg 

.0997 

.104 

15.08 

-.195 

.0138 

.1250 

.110 

15.37 

-.093 

.0168 

Mc=3.0 

12 

.0497 

.051 

15.15 

.248 

.0082 

.1006 

.061 

15  .74 

.459 

.0163 

5=10  deg 

- 

.1265 

.055 

15.39 

-.085 

.0208 

9 

.0374 

.050 

11  .47 

.908 

.0050 

A=0  deg 

.0742 

.033 

10.65 

-.213 

.0101 

.0929 

.045 

11  .23 

.214 

.0125 

Mc=0.4 

12 

.0374 

.969 

10.45 

-.112 

.0051 

.0749 

.986 

11.30 

.431 

.0104 

5=8  deg 

.0939 

.991 

11  .55 

.526 

.0131 

15 

.0373 

.889 

9  .86 

-.701 

.0040 

.0751 

.902 

10.38 

-.491 

.0086 

.0935 

.904 

10.46 

-.561 

.0108 

9 

.0374 

.026 

10.63 

-.054 

.0064 

A=0  deg 

.0746 

.033 

11.06 

-.369 

.0127 

.0934 

.042 

11  .61 

-.196 

.0157 

Mc=0  .4 

12 

.0376 

.979 

10.68 

-.008 

.0065 

.0749 

.000 

12.00 

.565 

.0131 

5=10  deg 

.0935 

.998 

11  .87 

.062 

.0163 

9 

.0502 

.138 

15 .10 

-.185 

.0045 

A =30  deg 

.0756 

.144 

15  .29 

-.148 

.0065 

.1003 

.149 

15  .44 

-.147 

.0083 

Mc=0.3 

.1259 

.145 

15  .32 

-.421 

.0108 

12 

.0498 

.065 

15.18 

-.103 

.0064 

5=8  deg 

.0656 

.074 

15.59 

.212 

.0082 

.0750 

.071 

15  .45 

.015 

.0095 

.0996 

.079 

15  .81 

.227 

.0122 

.1244 

.074 

15  .59 

-.142 

.0155 

.0503 

.060 

14  .94 

-.346 

.0065 

.1003 

.077 

15  .72 

.133 

.0124 

15 

.0503 

.013 

15  .65 

.364 

.0070 

.0758 

.009 

15  .45 

.011 

.0106 

.1000 

.018 

15.90 

.315 

.0139 

.1246 

.019 

15.95 

.216 

.0173 
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TABLE  4  (Concluded) 


A=30  deg 


M  =0.3 
c 


3=10  deg 


A=30  deg 


Mc=0.4 


3=8  deg 


A=30  deg 


Mc=0  .4 


3=10  deg 


aM 

( deg) 

T 

a 

e 

A 

(  deg ) 

( deg) 

9 

.0498 

.102 

14  .98 

-.111 

.0070 

.0750 

.110 

15.37 

-.013 

.0101 

.0994 

.114 

15.57 

-.095 

.0131 

.1241 

.115 

15.61 

-.341 

.0163 

12 

.0497 

.049 

15  .03 

-.060 

.0083 

.0754 

.050 

15.09 

-.297 

.0125 

.0995 

.060 

15.68 

.014 

.0162 

.1241 

.066 

16.03 

.079 

.0198 

15 

.0499 

.005 

15.31 

.218 

.0087 

.0754 

.008 

15.50 

.113 

.0131 

.0995 

.010 

15.63 

-.036 

.0173 

.1239 

.018 

16 .13 

.181 

.0215 

.0983 

.016 

16  .00 

.348 

.0171 

9 

.0377 

.051 

11  .52 

-.269 

.0050 

.0502 

.055 

11  .71 

-.223 

.0066 

.0754 

.050 

11  .47 

-.752 

.0100 

.0945 

.069 

12.36 

.033 

.0107 

.1002 

.068 

12.32 

-.187 

.0127 

12 

.0376 

.003 

12.15 

.362 

.0052 

.0496 

.010 

12.50 

.574 

.0069 

.0745 

.014 

12.70 

.488 

.0104 

.0933 

.016 

12.80 

.372 

.0130 

.0992 

.015 

12.75 

.254 

.0138 

.0377 

.998 

11  .90 

.111 

.0053 

.0753 

.002 

12.10 

-.121 

.0105 

15 

.0376 

.929 

11.55 

-.238 

.0047 

.0499 

.934 

11  .78 

-.149 

.0064 

.0751 

.941 

12.10 

-.119 

.0098 

.0944 

.954 

12  .72 

.280 

.0126 

.0997 

.943 

12.20 

-.301 

.0130 

9 

.0371 

.035 

11 .18 

-.107 

.0063 

.0495 

.042 

11  .61 

.096 

.0083 

.0743 

.045 

11  .79 

-.176 

.0125 

.0929 

.050 

12.09 

-.216 

.0154 

.0990 

.052 

12.21 

-.207 

.0164 

12 

.0375 

.989 

11.31 

.016 

.0065 

.0497 

.992 

11  .50 

-.017 

.0087 

.0746 

.003 

12.19 

.218 

.0130 

.0935 

.000 

12  .00 

-.317 

.0163 

.0992 

.010 

12.63 

.209 

.0173 

15 

.0374 

.935 

11  .03 

-.263 

.0060 

.0494 

.949 

11  .85 

.338 

.0082 

.0740 

.953 

12.09 

.129 

.0124 

.0934 

.957 

12  .33 

.015 

.0157 

.0991 

.963 

12.70 

.281 

.0168 

%o=14-51 


M  =1  1  .576' 
c 


“vo'11'36 


Mc  =  ll  .487° 


—  i 


Fig.  1  Airfoil  Cross  Section  Showing  Chordwise  Measum 


CHORDWISE  COMPONENT  Or 
FREESTREAM  VELOCITY 


Co  •  CN  SIN  o  -  CC  cos  a 
Ct  •  CN  COS  a  ♦  Cc  SIN  a 


Fig.  2  Airfoil  Force  Schematic 


UNSTEADY  PRESSURE  COEFFICIENT  Cp  (*)  UNSTEADY  PRESSURE  COEFFICIENT 
(SUCTION  SURFACE)  (SUCTION  SURFACE) 


FIRST  HARMONIC  COMPONENT  FIRST  HARMONIC  COMPONENT 

PHASE  ANGLE  OF  Cp  (*)  AMPLITUDE  OF  Cp  (x) 


aM  =  9  deg 


0  0.2  0.4  0.6  0.8  1.0 

DIMENSIONLESS  CHORDWISE  POSITION,  x 


Fig.  4  Example  of  Chordwiaa  Praasura  Distribution  Comparison  —  Discrete  Points 
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UNSTEADY  PRESSURE  COEFFICIENT,  C 


A  =  0  deg 


-  UPPER  SURFACE 

- - LOWER  SURFACE 


A  =  30  deg 


DIMENSIONLESS  CHORDWISE  POSITION,  x 


Fig.  5  Example  of  Chordwise  Pressure  Distribution  Format  for  Comparison 

of  Adjacent  Columns 
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-  9  deg  orM  -  9  deg 

nr  -  16  S  deq  '  «=  170  deq 

(PEAK  ANGLE 
OF  ATTACK) 
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DIMENSIONLESS  CHORDWISE  POSITION, 


THREE  HARMONIC  COMPONENT  RECONSTRUCTION  OF 
PRESSURE  COEFFICIENT,  CD  (*) 
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Pressure  Carpet  Plot,  Leading  Edge  Region 
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Fig.  9  Pressure  Contour  Map,  Leading  Edge  Region 
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Fig.  10  Pressure  Carpel  Plot,  Entire  Chord 
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DINENOXONLES8  TIME,  FRACTION  OF  PERIOD, 
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Fig.  13  Illustration  of  Segmented  Wave  Speed  Near  Leading  Edge 
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Fig.  14  Illustration  of  Segmented  Wave  Speed  Over  Entire  Blade 
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iWaWb 


0  0  02  0  04  0.06  0  06  0  10  0.12  0.14  0  0  02  0  04  0.06  0.08  0  10  0  12 

REDUCED  FREQUENCY.  kc 


Fig.  16  Effect  of  Pitching  Amplitude  on  Overall  Wave  Speed 


OVERALL  AVERAGED 


0  3C  O  A  □ 


REDUCED  FREQUENCY,  kr 


Pig.  17  Effect  of  Mach  Number  on  Overall  Wave  Speed 
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Fig.  18  Effect  of  Sweep  Angle  on  Overall  Wave  Speed 


Fig.  19  Effect  of  Sweep  Angle  on  Weve  Speed  at  Several  Chordwlse  Locations 
for  MP «0.30  and  a«8' 


it 


'  =  3C‘#  A 


0  0.02  0.04  0.06  0.08  0.10  0  12  0  0.02  0  04  0.06  0  08  0  10  0  12 

REDUCED  FREQUENCY.  kr 


0  Effect  of  Swoop  Anglo  on  Wovo  Spood  ot  Sovorol  Chordwito  Locations 
for  Mr.  >0.40  and  a  *8° 


WAVE  SPEED,  V 


DIMENSIONLESS  CHORDWISE  POSITION,  x 


Fig.  21  Chordwlse  Variation  of  Local  Wave  Spaad  for  aM  « is  dag,  Mc  -0.40, 
and  o«8  dag 
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Fig.  22  Effect  of  Sweep  Angie  on  Vortex  Inception  Angle  Based  on  initial  Wave  Propagation 
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Linear  Least  Squares  Scatter  Versus  Pitch  Rate 


FROM  GORMONT  FOR  V230I0-1  58  AIRFOIL  AT  M  O  4 


Correlation  of  Relative  Vortex  Inception  Angle  with  V^T 


RELATIVE  VORTEX  INCEPTION  ANGLE 


VORTEX  INCEPTION  PITCH  RATE,  Av 


Fig.  27  Correlation  of  Relative  Vortex  Inception  Angle  with  Av 
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(a)  =  0  deg  (b)  <vM  =  9  deg 


Fig.  28  Mach  Number  and  Frequency  Similarity  Analysis  of  First  Harmonic  Component  Response 
Distribution  for  30  deg  Swept  Airfoil  (Pitching  Amplitude,  n  =8  deg) 


(a)  nrM  =  0  deg  (b)  <irM  -  9  deg 
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eep  Angle  Similarity  Analysis  of  First  Harmonic  Component  Response  Distribution 


UNSTEADY  PRESSURE  COEFFICIENT,  Cp  (*) 
(SUCTION  SURFACE) 


=  A  Hz  (kc  =  0  050‘ 

=  t  Hi  (k  =  0  100. 
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(UPSTROKE 
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o  =  8  deg 


4  0 

DIMENSIONLESS  CHORDWISE  POSITION,  x 


Effect  of  Frequency  on  Complete  Unsteady  Chordwfse  Pressure  Distribution  at  the 
Lowest,  Mean,  and  Highest  Angle  of  Attack  Points  of  Motion  Cycle 
(c*m  *  0  deg,  Mc  *  0.30, 5  ■  8  deg,  A  «  0  deg) 


UNSTEADY  PRESSURE  COEFFICIENT,  Cp  (x) 
(SUCTION  SURFACE) 


Fig.  31  Effect  of  Mach  Number  on  Complete  Unsteady  Chordwise  Pressure  Distribution  at 
the  Lowest,  Mean  end  Highest  Angle  of  Attack  Points  of  Motion  Cycle 
(a^  *0  deg,  f  "8  Hz,  a  *8  deg,  A  *0  deg) 
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UNSTEADY  PRESSURE  COEFFICIENT.  Cp  (*) 
(SUCTION  SURFACE) 


1 


Fig.  32  Effect  o!  Sweep  on  Complete  Uneteedy  Chordwise  Pressure  Distribution  et  the 
Lowest,  Mean,  end  Highest  Angle  of  Attack  Points  of  Motion  Cycle 
■  0  deg,  kc  *  0.100  (f  *  8  Hz,  Mc  *  0.30),  and  a  «  8  deg 
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ALPHA 
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0  045 


UNSTEADY  PRESSURE  COEFFICIENT,  C 


A  =  0  deg 


UPPER  SURFACE 
LOWER  SURFACE 


A  =  30  aeg 


Fig.  35  Chordwlse  Pressure  Results  (or  the  Unswept  end  Swept  Configure  (Ions  et  aM  «fl  deg, 

kc  *0.100  (f  *8  Hz,  Mc  »0.30),  end  5  *8  deg 
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DIMENSIONLESS  CHORDWISE  POSITION.  X 

Fig.  37  Upp«r  Surface  Chordwiae  Prastura  Behavior  of  Swept  Airfoil  at 
a  M  =  0  deg,  kc  *  0.075,  end  a  *  8  deg 


GKfl 


Pig.  38  Upper  Surface  Chordwise  Pressure  Behavior  of  Swept  Airfoil  at  ajyj  =9  deg. 

kc  =  0.075,  and  5=8  deg 
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Fig.  39  Effect  of  Frequency  on  Upper  Surface  Response  at  Mean,  Maximum,  and  Minimum 
Points  of  Motion  Cycle  (q^  *  9  deg,Mc  «  0.30,  5  =  8  deg) 
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Effect  of  Mach  Number  on  Upper  Surface  Response  at  Mean,  Maximum,  and 
Minimum  Points  of  Motion  Cycle  (a^  =  9  deg,  f  =  8  Hz,  5*8  deg) 


Fig.  41  Effect  of  Frequency  on  Upper  Surface  Leading  Edge  Region  Time  Histories  ( M_  =  0.30.  <v  =  8  deg.  A  =  0  deg) 


ALPHA 


82 


n  Upper  Surface  Leading  Edge  Region  Time  Histories 


THREE  HARMONIC  COMPONENT  RECONSTRUCTION  OF  PRESSURE  RESPONSE.  C 


DIMENSIONLESS  CHORDWISE  POSITION,  X 


Fig.  47  Effect  of  on  Swoop  Dependence  of  Upper  Surface  Roeponeo  ot  a 
Upstroke  (f  *  8  Hz,  a  ■  6  deg) 
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Fig.  S3  Steady-State  Lift  Coefficient  for  NACA  0012  Airfoil  at  Mc  =  0.40 
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